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Abstract 

We study an exactly solvable quantum field theory (QFT) model describing inter- 
acting fermions in 2+1 dimensions. This model is motivated by physical arguments 
suggesting that it provides an effective description of spinless fermions on a square 
lattice with local hopping and density-density interactions if, close to half filling, the 
system develops a partial energy gap. The necessary regularization of the QFT model 
is based on this proposed relation to lattice fermions. We use bosonization methods to 
diagonalize the Hamiltonian and to compute all correlation functions. We also discuss 
how, after appropriate multiplicative renormalizations, all short- and long distance 
cutoffs can be removed. In particular, we prove that the renormalized two-point func- 
tions have algebraic decay with non-trivial exponents depending on the interaction 
strengths, which is a hallmark of Luttinger-liquid behavior. 



1 Introduction 

In this paper, we study a QFT modeQ describing interacting fermions in two spatial di- 
mensions (2D). We will refer to this as the Mattis model since it is similar to a model 
first proposed by Mattis [1] as an exactly solvable model potentially relevant for 2D lattice 
fermion systemsll In recent work, we presented specific physical arguments that the Mattis 
model provides an effective description of spinless fermions on a square lattice with local 
density- density interactions if, close to half filling, parts of the underlying Fermi surface 
have no low-energy excitations [21 El H]. This proposed relationship gives natural short- 
and long distance cutoffs for the Mattis model, and it fixes the parameters of the Mattis 
model in terms of the lattice model parameters [31 H]. The aim of this paper is to present 
a detailed and mathematically precise construction and solution of the Mattis model. 



*Elcctronic address: jodw02@kth.se 
t Elect ronic address: langmann@kth.se 

By this we mean a quantum model with infinitely many degrees of freedom. 
2 The precise relation of our model to Mattis' is explained in Section [SJ Remark 8. 
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1.1 Fermi arcs 



The derivation of the Mattis model from a lattice fermion system in [51 H] is not rigorous in 
the mathematical sense but relies partly on physical arguments. The latter were inspired 
by arguments that have been applied successfully to other low-dimensional systems; for 
example, the Luttinger model [5j [6j [7J |8] can be derived in a similar manner from a system 
of spinless fermions on a one- dimensional lattice. A mathematically inclined reader is free 
to ignore the details of this physics background and take the Mattis model at face value: 
as an, in our opinion, interesting exactly solvable model for interacting fermions in 2+1 
dimensions. We therefore give in the following only a short, intuitive description of the 
Mattis model and its relation to lattice fermions. Readers interested in further details are 
referred to Appendix IA1 

The Mattis model describes a coupled system of four different flavors (r, s), r, s — ±, of 
spinless fermions. It consists of two sets of parallel ID Luttinger chains (corresponding to 
s — ±), with each chain consisting of left- and right movers (corresponding to r = ±). Its 
relation to the original lattice fermion model is best explained using Figure [TJ shown in the 
figure are four line segments in Fourier space that remain after a square Fermi surface has 
been truncated around its corners. We refer to these line segments as Fermi arcs. Each 
Fermi arc is associated with one-particle degrees of freedom identified by flavor indices (r, s) 
and momenta k. The truncated regions correspond to one-particle degrees of freedom that 
are assumed to be gapped by interactions and therefore are ignored. 



Figure 1: A truncated square Fermi surface in Fourier space consisting of four disjoint 
Fermi arcs labeled by (r, s) with r, s — ±. Indicated in the figure are also the definition 
of momenta k± near the Fermi arc (+, +) and the interactions of the (+, +) arc fermions: 
they interact with opposing (— , +) and neighboring (±, — ) arc fermions with couplings 
proportional to 71 and 72, respectively. The other interactions are fixed by invariance of 
the model under rotations by 7r/2. 
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The fermions thus are represented by field operators ^/v iS (k) in Fourier space, and the free 
part of the Mattis Hamiltonian has the form (see Section 12.11 for the precise mathematical 
definition) 

H o=Yl \ d ' k e ^( k ) ^l(k)^ s (k): (1.1) 

r,s=± J 

with dispersion relations linear in k s (orthogonal to each Fermi arc) and constant in k_ s 
(parallel to each Fermi arc), i.e. e r)S (k) = rvpk s with a model parameter vp > (referred to 
as Fermi velocity). The integration domain in (11.11) is unbounded in the fc s -direction, while 
\k-s\ < vr/a with 2ir/a the length of each Fermi arc. 

The interaction part of the Mattis Hamiltonian consists of density-density interactions 
between different Fermi arc fermions. There are two different coupling parameters: 71 for 
interactions involving fermions on opposite arcs, and 72 for fermions on nearby arcs; see 
Figure [TJ As seen in Section 12. 1[ the possible momentum exchange p in these interactions 
are restricted such that \p±\ < ir/a. Thus the length parameter a, which is proportional 
to the lattice constant in the original lattice fermion model, serves as a natural ultraviolet 
(UV) cutoff in the Mattis model. There is also a natural infrared (IR) cutoff provided by 
the linear size L of the original square lattice. 



1.2 Summary of results 

The Mattis model can be formally^ defined by the Hamiltonian 

r,s=± J 



(1-2) 

+ 9r,s,r',s> ■^l, s %s--^l',Ar',s''-) 

r',s'=± 

with d± = d/dx± and x± Cartesian coordinates of positions x. The ijjr} = ipr)(x) are 
fermion field operators obeying canonical anticommutator relations {ip r a (x), ijr r i s /(y)} = 
8 r y8 SjS i8(x — y), etc. As already mentioned, the model parameters include the Fermi 
velocity vp > 0, the coupling constants g r , s ,r',s' parameterized as 

g r ,s,r', s ' = a7ro F (6 SjS '5 r _ r /7i + 8 s ^ s >y 2 /2) (1-3) 

with dimensionless parameters 7^2, and the UV cutoff a. The colons indicate fermion 
normal ordering. 

Our first aim of this paper is to give a rigorous definition of the Mattis model; see 
Section 12.11 For that we consider the regularized version obtained in [3] (in the special 
case 71 = 72 > 0). This regularized Mattis model describes an infinite number of fermion 
degrees of freedom, but it comes with particular short- and long distance regularizations 
controlled by the UV- and IR cutoff parameters a and L, respectively. As we will see, this 
regularized Mattis model is well-defined if and only if 

|7i|<l, |7 2 |<|l + 7i|. (1.4) 



3 We suppress details of the UV regularization here. 



3 



We emphasize that this parameter restriction does not mean that the Mattis model is 
weakly coupled: as shown in [3], the scaling of the coupling (11. 3p naturally arises in a 
partial continuum limit of the above-mentioned lattice system, and 71 = 72 = 1/2, for 
example, does not correspond to weak coupling in the latter system. 

A key result (see Section 12. 2j) is to make mathematically precise the equivalence of the 
Mattis model to a non-interacting boson model formally defined by the Hamiltonian [3] 

H = f ^ / * x '■ ( (1 ~ 7l)n * + (1 + 7l) ^ $s)2 + 72 ^ $s) : ( L5 ) 

with 2D boson field operators U s = II s (x) and $ s = $ s (x) satisfying canonical commutator 
relations [IL,(x), $ s /(y)] = — i8 StS >5(x — y), etc. The boson operators are related to the 
fermions as 

d s $ s = v / vra(: ip+, a ip+ tS ■ + ■ ^-,s^-,s ■) ^ g v 

n s = -V^r5(:V'+ lS '0 + , s : - :ipl >s ip- )S -) 

at finite UV cutoff a. This relation, and the non-trivial scaling of the coupling constants 
gr,s,r',s' with a, make it clear that the UV cutoff in the Mattis model is essential. As we 
will see, the full continuum limit a — > + of the model, referred to as the UV limit in 
the following, makes sense only after non-trivial multiplicative renormalizations. Similarly, 
L — > 00 is referred to as the IR limit. 

Since the regularized Mattis Hamiltonian can be written as a system of free bosons, it 
can be diagonalized by analytical methods. The result is given in Theorem 13.11 We also 
compute the free energy (Result l3T2~j) . the density- density correlation functions (Result [33]) . 
and the fermion iV-point correlation functions (Result 13. 6p . We note that all these results 
are exact in the thermodynamic limit, and they depend on the UV cutoff a. Our results 
suggest that it is possible to remove the UV cutoff and, by appropriate renormalizations, 
obtain non-trivial results for all observables of the Mattis model. In the present paper, we 
give detailed proofs of this for the free energy and the fermion two-point functions (the 
latter at zero temperature); see Results I4.1H4.4I We find that these two-point functions 
in this limit have an algebraic decay with exponents that depend non-trivially on both 
coupling parameters 71^, which is one characteristic of Luttinger-liquid behavior [9] [TO] . 
We also discuss how our results can be used to substantiate the derivation of the Mattis 
model from the original lattice fermion system [3J; see Section Remark 5. 



1.3 Related work 

The present paper is the third in a series [3J H] in which we propose and develop a method 
to do reliable, non-perturbative computations in models of strongly interacting fermions 
on a 2D square lattice [2]. We have so far considered the simplest non-trivial case, the 
so-called 2D t-t'-V model, describing spinless fermions with local interactions. In [3], a 2D 
analogue of the Luttinger model was derived from the 2D t-t'-V model using a particular 
partial continuum limitE) This 2D Luttinger model is an extension of the Mattis model 

4 The interested reader can find more details in Appendix 1X1 
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by so-called antinodal fermions that cannot be bosonized. In [I] we showed that, within a 
mean field approximation, there is a non-trivial phase away from half-filling in which the 
antinodal fermions are gapped, and in this phase the 2D Luttinger model reduces to the 
Mattis model. One key ingredient of our approach is the notion of underlying Fermi surface 
arcs (not necessarily corresponding to a true Fermi surface) in the 2D t-t'-V model away 
from half filling. Previous theoretical work predicting four disconnected Fermi surface arcs 
in 2D interacting lattice fermion systems include renormalization group studies in jTH [12] . 

One of our main tenets is that progress in understanding these lattice models can be 
achieved if one clearly distinguishes between approximations justified by physical arguments 
and manipulations that are mathematically exact (rigorous); the present paper belongs to 
the latter category. We mention that the model derived in [3] corresponds to the special 
case 7 = 71 = 72 > 0; the generalization to different coupling parameters, allowing also for 
negative values, is natural from a mathematical point of view. We also note that, while an 
outline of how to diagonalize the Mattis Hamiltonian by bosonization was already given in 
[3], Section 6.2, the details and complete solution of the model are given here for the first 
time. 

We choose the name of the model studied in this paper to acknowledge a pioneering 
paper published by Mattis in 1987 in which he pointed out that a 2D interacting fermion 
Hamiltonian similar to ours can be mapped exactly to a non-interacting boson Hamiltonian 
PQ; the dispersion relations of this boson Hamiltonian were given by Mattis (see [1], Equa- 
tion (8)) but no details were provided. Mattis also argued that this model can arise from a 
tight-binding description of 2D lattice fermions at half-filling with a square Fermi surface 
and in the absence of nesting-type instabilities, but, as shown in [3], this latter proposal 
is doubtful (at half-filling, there are additional interaction terms that cannot be bosonized 
in a simple manner and which are likely to yield a gap, very different to what the Mattis 
model describes). 

Mattis' insight described above has received little attention up to now, and it was 
re-discovered in 1994 independently by Hlubina [T3j and Luther [T3]. Hlubina presented 
a model similar to Mattis' that he claimed can be mapped to a non-interacting boson 
Hamiltonian, and he diagonalized the latter. He also gave phenomenological arguments 
suggesting that this model describes a two-dimensional Luttinger liquid. However, from 
his discussion, it is not clear what the solvable model actually is (the model formulated 
in [13], Equation (1), is not exactly solvable but needs to be modified — the necessary 
changes are similar to ones discussed in detail in [3]). Luther studied spinfull fermions with 
a square Fermi surface and argued that the non-interacting part of the Hamiltonian can 
be simplified to one that can be bosonized exactly, and which is essentially the same as in 
Mattis' model. Luther diagonalized the bosonized Hamiltonian that one obtains by also 
including density-density interactions, initially restricting to interactions between opposite 
Fermi surface faces only [13] (this restriction corresponds to 72 = in our notation). This 
was later extended to also include interactions between adjacent faces in [T5] . 

Our solution of the Mattis model not only includes the diagonalization of the bosonized 
Hamiltonian, but also the computation of all fermion correlation functions using the boson- 
fermion correspondence formula, which is an operator version of a formula that first ap- 
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peared in [16]. In regards to this correspondence formula, we emphasize the role of the, 
at least in the condensed matter literature, sometimes overlooked Klein factors, i.e. the 
unitary operators that lower- or raise fermion excitation numbers when acting on states in 
the Fock space. For one- dimensional systems, the unitary Klein factors often combine to 
unity when computing physically interesting correlation functions, and they can then be 
safely ignored, as was done in the well-known work by Luther and Peschel [17]. However, 
the situation in higher dimension is different due to the profusion of flavor indices; previous 
work on computing correlation functions by Luther and collaborators [T8| IT5] neglected to 
incorporate Klein factors, and their results are therefore not directly applicable to the model 
studied in this paper. In particular, as will be shown, Klein factors cannot be ignored in the 
solution of the Mattis model. In the condensed matter literature, the importance of Klein 
factors has been pointed out in papers by Haldane [19j [9] and Heidenreich et al. [20] in the 
context of the Luttinger model; see also [21] for a review on constructive bosonization and 
Klein factors, and Remark IB. 51 in Appendix [Bj Our formulation and solution of the Mattis 
model is based on mathematical formulations of these results proved in [221 [231 [2H [2U 26J, 
for example. 

An important question is whether the Mattis model describes a 2D Luttinger liquid [10J. 
Different regularizations and treatments of Mattis-like models can give different results, as 
exemplified by Refs. [2Z] arid [28]: these works are on such models that differ in regulariza- 
tion details of the interaction, they use different methods, and they obtain different results. 
This suggests to us that it is important to be mathematically precise both in the definition 
of the model and of Luttinger-liquid behavior, as well as in the methods used to study 
them. 

There exist many other papers studying 2D fermion systems by mapping them to bosonic 
models; see e.g. [29] and references therein. However, these differ significantly in detail from 
our constructive bosonization methods. In particular, many of the various approaches, as 
indicated above, do not establish an operator correspondence between fermions and bosons, 
thus precluding the exact computation of fermion correlation functions. 

1.4 Notation and conventions 

We write "=" to emphasize that an equation is a definition. For z G C, the complex 
conjugate is written z, the real part is ^ft(z), the imaginary part is ^(z), and the argument 
(phase) is arg(z). We write a 2D vector as u = u + e + + w_e_, with e± • e± = 1 and 
e± e T = 0. If S is a set, x G yS and x G y + S means xy^ 1 G S and x — y G S, respectively. 
We identify 1 with the identity operator. Unless otherwise stated, and when there is no 
risk of confusion, a free variable or index in an identity can take any allowed value. For 
example, the relations in (12.21) below are true for all r, s, r', s' — ± and k, k' G A*, whereas 
in equation (12.31) . the identity holds true for all r, s — ±, but only those k G A* satisfying 
the given condition. 
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1.5 Plan of paper 

A detailed definition of the Mattis model, including all regularizations, is given in Sec- 
tion |2j This section also contains the mathematical results needed to bosonize the Mattis 
model, and a technical discussion on operator domains and the essential self-adjointness of 
the Mattis Hamiltonian (this latter part may be skipped without loss of continuity). In 
Section [21 we solve the Mattis model by diagonalizing the Hamiltonian and by computing 
all correlations functions. In Section HI we give the renormalized free energy and fermion 
two-point functions after all short- and long distance cutoffs have been removed. Section [5] 
contains final remarks. Technical details are deferred to Appendices lATlFl Appendix lAl gives 
further details on the relation of the Mattis model to lattice fermions. In Appendix [Bj we 
summarize well-known mathematical results on constructive bosonization in one dimen- 
sion, and we explain how these are used to prove the results in Section [2J Appendix O 
contains computational details of the solution of the Mattis model. Details on the removal 
of short- and long distance cutoffs in the Mattis model can be found in Appendix |D] In 
Appendix [Ej we discuss certain terms that were neglected in the derivation of the Mattis 
Hamiltonian from the lattice fermion Hamiltonian in [3]. Appendix |F] lists some index sets 
used throughout the paper. 

2 Prerequisites 

In this section, we define the regularized Mattis model in terms of fermion creation- and 
annihilation operators. We also summarize the mathematical results needed to bosonize 
the model. 

2.1 The Mattis model 

We introduce fermion field operators ipfi) (k) labeled by flavor indices r, s = ± and momenta 
k in the following sets^l 

, „ def [ , 2-7T / 1 \ 2 7T 7T 1 . 

A;=(k £T ( Z + -j (2.1) 

with parameters a > and L > such that L/a is an odd integer. These operators obey 
the canonical anticommutator relations 

{i iS (k),$ v (k')} = (^) V'V<W', {^, s (k),i v (k')} = (2.2) 

and are defined on a fermion Fock space T such that ^> rs (k) is the Hilbert space adjoint 

of ^ s (k). The representation of these operators on T is specified by the existence of a 
vacuum state (Dirac sea) Q G J 7 satisfying 

Vv iS (k)ft = $, a (-k)fl = 0, Vk such that rk s > 0, (2.3) 

J This and several other index sets are collected in Appendix IF1 for easy reference. 
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and (Q,Q) = 1, with (•, •) the inner product in J 7 . We define normal ordering of fermion 
bilinears O = ^| )S (k)^ r , jS ,(k') by 

:0: = O - (Q,OQ). (2.4) 
The same definition will be used for boson bilinears. 

The following characteristic (cutoff) function will be needed 

x(p) «/i if - */a < p± < -I- (2 5) 

otherwise 



We also introduce 

ki,k 2 SA* neZ 



4 S (P)= Yl (^") : ^( kl )^,s( k 2): ^5k 1+ p,k 2 +2 T n5-ie_ s (2-6) 



for p G (2tt/L)Z 2 and r, s = ±. We refer to J r , s (p) as (fermion) densities. We use the same 
name for the Fourier transform J r ,s( x ) of these operators (to be defined below). 

Definition 2.1. Let > and g r ,s,r',s' as i n HI -3j) (II .4p . Then the regularized Mattis 
Hamiltonian is given by 

H = H + H 1 (2.7) 

with the free part 

Ho = v F J2 E (t)^ : ^( k )VUk): (2-8) 

r,s=± keA| 

and the interaction part 

Hi = E E (r) &ns,r', S 'X(p)^r, S (-p)^r',8'(p)- (2.9) 
r,r',s,s'=± pe^LZ 2 



Our results in Section [3] prove that the regularized Mattis Hamiltonian is a self-adjoint 
operator bounded from below. 

Remark 2.2. The densities in ( 12.6P and the Hamiltonians in ( 12.7l) -( l2~9j) are all unbounded 
operators, and a complete definition of the Mattis model thus require both a specification 
of their common domain of definition V, and a proof that Hq and H are essentially self- 
adjoint on V [30j [31]. This can be done using general results on operator algebras arising 
in quasi-free boson- and fermion quantum field theories [32]; the interested reader can find 
the details in Section 12.31 An informal description of this goes as follows: In Lemma 12.7( b) 
below, we give a complete orthonormal set of exact eigenstates of H , and this proves 
essential self-adjointness of -f^o o n the vector space T> of finite linear combinations of these 
eigenstates. Moreover, it is easy to see that D is a common, invariant domain for the 
operators in (I2.6p - fl2.9l) . and for other operators to be introduced below. For this reason, 
all identities involving commutators of unbounded operators that we write down are well- 
defined as identities on V. Strictly speaking, one should distinguish between an unbounded 
operator A defined on V and its closure A (see e.g. [30], Section VIII. 1), but, as can be 
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inferred from Section 12. 3[ it is safe to abuse notation and use the same notation for both 
operators. Finally, one can prove directly that the Mattis Hamiltonian H is essentially 
self-adjoint on T> (Corollary 12. lip , but, as already mentioned, this fact is also implied by 
our exact solution of the Mattis model in Section [3] (in which a complete orthonormal set 
of exact eigenstates of H, all contained in D, is found). 

We can write the regularized Mattis model in position space, thereby making precise the 
formal Hamiltonian given in the introduction. Our conventions for the Fourier transform 
of fermion field operators are 

<M x ) = h £ (x) 2 ^ s(k) eik ' x (2 - 10) 

keA* 

with x in the following "position space" corresponding to (12. ip . 



A, = { x : x s E R, x_ s G aX, ~ < x ± < ^} . (2.1 I 



which is different for s = + and s = —. Note that this is not a subset of the lattice in the 
original fermion model [3]: one component in x 6 A s takes on continuous values, the other 
discrete ones, and the "lattice constant" a differs from the original one. 

The operators in (I2.10p satisfy the anticommutation relations 

{W, s (xUl v (y)} = <5 r ,A,A( x -y) (2-i2) 

with the partially regularized Dirac delta function 

^(x) d = f 5(x s )^_ s ,o. (2.13) 
a 



Likewise, the Fourier transform of the densities in (12. 6 p are defined aa, 

fr,s(x; e) 

!\ 2 f ^„ip.x„- e |pJ/2 n , (2-14) 



J r;S (x) = lim J riS (x;e) 



J r , s (x; e) d = f (l) lM e ip ' X e" e|Ps|/2 (e > 0) 
P eA* 

with 

A; = (pG> 2 : -Z<P-.<l). (2.15) 
I L a a J 

One finds that 

J P , a (x) = :^, a (x)V M (x):, (2.16) 
i.e. the fermion densities are local operators. We introduce the following notation 

f r L / 2 ^ 

d 2 x = / dx s J2 ~ a (2- 17 ) 



'Only J rjS (x) is needed in this section, while J r , s (x; e) for e > is used in Section 
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with a ID lattice 

Aid = |x G aZ : -~ < x < ^ j (2.18) 

for the discrete component in (12.111) . The regularized Mattis Hamiltonian can then be 
written in "position space" as (cf. (II. 2p ) 

r s s 

(2.19) 

+ flW'.s' / d 2 y J ns (x)5(x -y)J r / jS /(y) ) 
with the fully regularized Dirac delta function 

^( x )= E (i) 2 *(p) e " ip ' x ( 2 - 2 °) 

Remark 2.3. For future reference, we note that ^v iS (x) in (I2.10p and J r ,s(x) in (12. 14|) are 
well-defined for all x G M 2 , although, to invert ( 12 . 1 Of) . only ip r>s (pc) for x G A s is needed: 
Vv,s(k) = J s d 2 x , ?/' r . jS (x) e~ lkx /(27r), and similarly for J r>s (x). Later, when considering cor- 
relation functions in Section |3] and HJ we find it convenient to consider the variable x as 
continuous. 

Remark 2.4. For simplicity in notation, we use in this paper a UV regularization that differs 
somewhat from the one obtained in [3]. We expect that details of the UV regularization 
can be changed without significantly changing physically relevant results. Moreover, the 
generalization of our results to the more complicated UV regularization in [3] is, in principle, 
straightforward. In any case, both regularizations coincide in the special case k — 1/2 (see 
Appendix lAl for the definition of k). 



2.2 Bosonization 

In this section, we collect the mathematical results needed to bosonize the Mattis model. 
Some additional details can be found in Appendix [0 

Proposition 2.5. (a) The fermion densities are well-defined operators on T obeying the 
commutation relations 

Jr,s{p), Jr',s'(P ) 

neZ 

Moreover, 

X, s (p) t = IA-P) (2-22) 

To be precise: These operators and the identities below are well-defined on a common, dense, invariant 
domain T> to be defined in Section |2~51 see also Remark 12.21 



rS rr '8 ss > f S (—) 'Wp', 2 7 rng- 1 e_ s ■ (2-21) 

a, \zir/ * — ' 



10 



and 

J rs (p)Q = 0, Vp such that rp s > 0. (2.23) 
(b) The following operator identity holds true, 

E (t)^ : ^( k )^, s (k): = arr £ (I)' : J r ,.(-p) J r ,.(p) : (2.24) 
un'i/i 6ot/i szrfes defining self-adjoint operators on T (the colons indicate normal ordering as 



{Proof: See Appendix iBl) 

Remark 2.6. It follows from equations (I2.14p and (12.21)) that 



1 d ~ 

J r , a (x), Jr', S '(y)] = ^rAy^— ^— 5 s (x - y). (2.25) 

27rai ax s 



This and ( 12 .16)) imply the identities in (II. 6p . 

From Proposition 12.5( a) follows that the rescaled densities 

-i./gf4 s (p) if Ps >0 



and 6j(p) = & s (p)*, for p in 

K = {p G A; : p s ^ 0} , (2.27) 

obey the defining relations of boson creation- and annihilation operators, 

[& s (p),^(p')] = ^(p)A'(p')] = 0> ^(P)^ = 0. (2.28) 

Furthermore, 12.5( b) implies that (12.81) can be written as 



H = aTrv F E E (I)' : 4*(-P)4-(p) : ( 2 - 29 ) 

r,s=± p6 Aj 

and thus the regularized Mattis Hamiltonian can be represented solely in terms of densities. 

It is important to note that not all Fourier modes of the fermion densities can be 
expressed as bosons: there are also the zero modes 

4.(P-.) = \/J- MP) (2-30) 
V ^vr p s =o 

with p_s in 



AId = <|pG^Z : --<p<-\. (2.31! 
L a a 
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(2.34) 



which require a different treatment. We find it convenient to consider the following Fourier 
transformed zero-mode operators 

N r>s (x) = l^rM eipx i x e Aid- (2.32) 

It follows from that 

N rjS (x)tt = 0. (2.33) 

To fully formulate the Mattis model in bosonized language, we also need to introduce the 
Klein factors R r , s (x) conjugate to the zero modes in (I2.32p . These, together with the boson 
operators, span the full fermion Fock space J 7 : 

Lemma 2.7. (a) There exist unitary operators R r , s (x) on the fermion Fock space J 7 com- 
muting with all boson operators in 112. 26]) and satisfying the commutation relations 

[N r)8 (x),R r > j8 >(x')] = r8 rjr >8 SjS >5 XiX >Rr tS (x), 

(b) Let Q be the set of all pairs (n, v) with 

n = R(p)} s=±jPe A s * > ,/ = K.(«)} r , a= ± i<B6AlD ( 2 - 35 ) 
and integers u r>s (x) and n s (p) > such that 

E E u r,s( x ) 2 < °°j E E I^I^(p) < 00 • ( 2 ' 36 ) 

r,s=± :reAi D s=± pg A* 

JTien i/ie states 

= (n n ^S) ( n n ^w^ w ) n - (2 - 37) 

wif/j (n, i/) G Q, provide a complete orthonormal basis in J 7 . 

(c) The states r] n u are common eigenstates of N r>s (x) and 6j(p)6 s (p) with eigenvalues 
u r)3 (x) andn s (p), respectively. 

(Proof: See Appendix iBl) 

The above definitions allow to rewrite H in ( I2.29P as 

Ho=v F j2J2 ipMipMp) + ^E E N ^ 2 ( 2 - 38 ) 

s=± p( z\* r,8=± xeAm 

where we used the Parseval-type relation 

peA* 1D " xeAm 

This and Lemma 12.7( c) imply that the states in ( 12. 37ft are exact eigenstates of Ho with 
eigenvalues vp J2 S \p s \n s (p) + nvp s x f r>s (x) 2 j L, and the latter are finite due to ( I2.36p . 

To compute fermion Green's functions, we need the following result expressing the 
fermion operators in terms of the bosons: 
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Proposition 2.8. For s = ±, x G A S; and e > 0, the operators 

\Jlixae 

xex J r ± V f^V-X.sfple^e-^l/ 2 



eXT) ( r ±T (-Y-J (P) - - ^ A (2 ' 40) 

ex p^ 27r Ps JrAP) 

P eA* 

are proportional to unitary operators on the fermion Fock space Moreover, 

Vv jS (k) = lim ±- [ d 2 x V r , s (x; e) e" ikx (2.41 

e-5-0+ Z,7Tj s 

using the notation in {2. i7| j|^| 



(Proof: See Appendix iBl) 

Note that, by taking the Hilbert space adjoint, (I2.4ip implies a corresponding result for 

The operator V riS (x;e) can be regarded as a regularized version of the fermion field 
operator t/v iS (x). To be precise: ^/v )S (x) is an operator- valued distribution, and products 
of these need to be interpreted with care. However, for e > 0, V rjS (x;e) is a well-defined 
bounded operator, and it converges to Vv,s( x ) m the limit e — > + . This also suggests 
to write ?/v,s(x;e) instead of V r)S (x;e). In Section [373| we compute expectation values of 
products of fermion field operators as the limit e — > + of the corresponding expectation 
values of products of these regularized operators. This circumvents potential difficulties 
with the distributional nature of the fermion fields. 

We find it convenient to introduce the notation 

Q r ,.(x) =4=(V + , s (x) + rN-,.(x)) 

Y V \ (2.42) 

QrAp)=^(N + Ap)+rN^ s (p)) 

In our computations in the next section we also need the following orthogonality result 
implied by Lemma 12. 7\ 

(Vn,u, J] RrAx) mrAX) Vn,v) = J] <W (*),(> (2-43) 
r,s,x r,s,x 

independent of (n, u), for all integers m rjS ( x ) such that J2 r s x m r,s( x ) 2 < 00 an d all states 
Vn,u in (I2.37P ; the products in ( I2.43P are over all r, s — ± and x G A 1D . Note that the 
product on the left-hand side in (I2.43P is only well-defined with an ordering prescription 
(since operators R r;S (x) and R r i^(x') anticommute for (r, s,x) ^ (r',s',x')). However, the 
result in (I2.43[) holds true for any ordering, and we therefore do not specify the latter. 



8 To be precise: v27n5eV r , s (x; e) is a unitary operator. 
9 The limit is in the strong operator topology; see e.g. [50] , 
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2.3 Unboundedness of operators 

In general, sums or products of unbounded Hilbert space operators need not be densely 
defined, and a symmetric unbounded operator need not define a self-adjoint operator. How- 
ever, as elaborated in this section, such difficulties do not arise if a set of unbounded opera- 
tors belongs to what is called a Op* -algebra, and this is the case for the unbounded operators 
introduced above. We also prove that the Mattis Hamiltonian is essentially self-adjoint on 
V. 

We recall the definition of a Op*-algebra (see |32j, Section 11(c) and references therein). 
Let D be a dense subset of some Hilbert space (in our case J 7 ). Then C + (T>) is the *- 
algebra of linear operators with T> as common, invariant domain of definition, and with the 
involution + defined as A + = A*\-p, i.e. A + is such that 

(771, Ari 2 ) = (A + rn, rj 2 ) V771, rj 2 E V (2.44) 

(in this section, we also use the symbol * to denote the Hilbert space adjoint). It follows 
from this definition that all elements A e C+(T>) are closable [30]. An Op*-algebra AonV 
is a *-subalgebra of £ + (T>) containing the identity. It is called standard if each symmetric 
element in A is essentially self-adjoint. 

We define T> as the space of all finite linear combinations of vectors in f 12 . 3 7[) . By 
definition, this set is dense in J 7 and a domain of essential self-adjointness of Hq. The 
densities J rjS (p), the free fermion Hamiltonian H Q , and the regularized Mattis Hamiltonian 
H, are all unbounded operators defined on the common, invariant domain T>. It follows 
from the above definition that these operators and the Klein factors R rtS (x), generate a 
Op*-algebra containing the Mattis Hamiltonian H. This makes mathematically precise the 
identities involving unbounded operators given in the previous section; they are well-defined 
as Op*-algebra relations on the domain D if we interpret t as involution + (rather than as 
Hilbert space adjoint). A non-trivial question is if this Op*-algebra is standard. However, 
this need not be addressed here in full generality. Instead, we restrict ourselves for the 
remainder of this section to a particular subspace of the Op*-algebra, which contains the 
Mattis Hamiltonian, and for which the proof that all symmetric elements are essentially 
self-adjoint becomes simpler. 

For < E < 00, let P E = 6{E - H ) with 6(x) = 1 for x > and for x < 0, i.e. P E is 
the projection to the Hilbert subspace generated by all eigenstates of H with corresponding 
eigenvalues less or equal to E. We define the set O u (Ho) of all linear operators A on T 
such that, for all E < 00, AP E is bounded and 

\\AP E \\ < a{E + E ), AP E = P E+A AP E (2.45) 

(|| ■ || is the operator norm) for some finite and non- negative integers a, Eq and A (all 
independent of E). Note that O w (Hq) is a vector space but not an algebra. 

Lemma 2.9. Symmetric operators A e 0"(Hq) are essentially self-adjoint on T>. 
Proof: For all rj G T>, there exists E < 00 such that rj = P E rj, and thus, using (I2.45p . 

n-l 

\\A n r)\\ = \ \AP E+{n _ 1)A - ■ ■ AP E+A AP EV \ \ < Y[a{E + jA + E )\\ v \\ (2.46) 

3=0 
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implying J2 n =i \ \A n r]\\t n / n\ < oo for < t < l/(aA) by the ratio test. Thus V is a set of 
analytical vectors for any A G O w (Hq), which gives the result by Nelson's analytic vector 
theorem (see e.g. [31], Section X.6). □ 

Remark 2.10. More generally, A G O u {Hq) implies A + = A*; see [32], Lemma 1. 
Corollary 2.11. The Mattis Hamiltonian H is essentially self-adjoint on T>. 

Proof: We only need to prove that H G C^(ifo)- The commutator relations 

[H , J r ,s(P)] = - rV FPsJr,s(P) ( 2 - 47 ) 

imply the second condition in (I2.45[) for A = H, with A = l-nvpja. 

To prove the first condition in (12.451) we consider A = J r ,s(—p) X;s(p) ■ By definition of 
normal ordering, A = C\+ : A : with C\ = rp s 8(rp s )L 2 j (2ira) . Thus, for arbitrary rj G D 
and E < oo such that rj = PeVi 

(rj, Arj) = (V, (ci+ : A :)t)) < ( V , ( Cl + c 2 H )ri) < (c x + c 2 E) (2.48) 

with C2 = L 2 /{a/KVp). This and the Cauchy-Schwarz inequality imply the first condition in 
(I2.45P for all operators A = J r - yS {— p)J r ',s'{p)- Since the Mattis Hamiltonian is a sum of H 
and a finite linear combination of such operators A, the triangle inequality and | \HqPe\ \ < E 
imply the result. □ 



3 Solution of the Mattis model 

In this section, we show that the free energy and any correlation function of the Mattis model 
can be computed by analytical methods. Computational details are given in Appendix O 

3.1 Diagonalization of Hamiltonian 

Proposition 12.5( a) implies that, for s = ± and p G A*, the following linear combination of 
fermion densities 

y^ sV ; ( Ps ^o) (3.i) 

n s (p)=y£(-J + , s (p) + J-,,(p)) 

obey the canonical commutator relations of 2D neutral bosons, i.e. 

[^(p),n;,( P ')] = i^(^) 2 w> [n s ( P ),n;,(p')] = [*.(p),&M = o, (3 2) 
n+(p) = n.(-p), £t(p) = $ (_ p ). 
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Moreover, Definition 12.11 and Proposition 12.5( b) imply that the Mattis Hamiltonian can be 
written in terms of these aa 10 l 



H = y E E (t) 2 : ( (1 - lix)fl ' A + ™)rf* f A 

s=± P eA* (3.3) 

with 

2 s (p) = -72iP«x(p)Q+,-«(p«)^-.,o, Sj(p) = S a (-p) (3.4) 
the zero modes that couple to the bosons, and 

f q ) = t(t) 2 (E E (i + ni)Q r A-p)QrAp) + ^2Q + A0)Q + A0)) (3.5) 



r,s=± 



PGA* r 



the terms containing only zero modes; we suppress common arguments p in ( I3.3p . Here 
and below, terms labeled by U B" and "Q" generically correspond to boson- and zero-mode 
contributions, respectively. 

The zero- mode terms in ( 13. 4ft and (13.5j) commute with all boson operators Il s and $ s , 
and they can therefore be treated like C-numbers. It is straightforward to find a Bogoliubov 
transformation that diagonalizes the Hamiltonian in (13.31) into a sum of decoupled harmonic 
oscillators and zero-mode terms. In the following we summarize the main results from this 
diagonalization; see Appendix O for details. 

Theorem 3.1. There exists a unitary operator IA on the fermion Fock space J 7 diagonalizing 
the regularized Mattis Hamiltonian as follows: 

H d ^U^HU = H B + H Q + £ 

= E E w.(p)6!(p)6.(p) l H Q = H Q ({QrA*)}) (3 - 6) 

«=± pGA* 

with the boson operators b s (p) in $2.26)) - $2.23) and 



^, + {p) = ^ ^VKlPl 2 ±>/|p| 4 -A(2p + p_)2) z/ 7 2X(pKp- ^ (3 ?) 

.ufa/C 1 -7ix(p))b±l l f 12X(P)P+P- = o 



A - 1 -(i^- 1 ) • 'V='vVi-'r (3.8) 



£/ie boson dispersion relations, 

H Q ({Qr,s(x)}) = ^J2 E ((l+7i)^ + , s W 2 + (l-7i)Q-, s (^) 2 ) 



± xeAi D 



- + 7i 

=± xsAid xsAid xeAio 



(3.9) 



10 Note the slight abuse of notation: the summation over momenta has been extended such that it super- 
ficially includes boson operators that are not defined. However, all these operators come with vanishing 
prefactors. The same remark applies to various expressions in Appendix [Q 
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the zero-mode contribution with the charges Q r , s {x) in \2.30\) - &2.32\) and ^Tj2 ), and 



the groundstate energy. 



(Proof: See Appendix O) 



(3.10) 



s -± pGA| 



The explicit expression for the unitary operator U is given in (1C33j) . (1C.39j) . and (1C.40]) . 
Note that the sum in (I3.10P has only a finite number of non-zero terms. Furthermore, the 
conditions in (ll.4p are necessary and sufficient for the dispersion relation in (I3.7p - (l3.8p to 
be well-defined (the conditions imply that < A < 1). 

When 72X(p)p+P- 7^ 0, the boson dispersion relations w±(p) in (13.71) are given by vp\p\ 
times dimensionless functions depending only on the polar angle 9 = arctan(p_/p + ). The 
angular dependence is determined by the constant A in (13.81) . and it is qualitatively similar 
for all values of A. For future reference, we note that both functions U}±(p)/(vp\p\) are 
invariant under 6 — > 9 + tt/2 and 9 — > —9, and oj-(p) / (vf\p\) vanishes like const. \9\ as 
9 -> 0, for all |p| > 0; see Figure EJ 



Q_ 

LL 
> 



0.25 



V) 




Figure 2: Solid line: the angular dependence of the boson dispersion relations w+(p) (top 
curve) and w-(p) (bottom curve) in ( 13. 7p for fixed |p| > and 71 = 72 = 0.9. Dotted line: 
same angular dependence using instead the effective dispersion u)±(p) defined in (14. 2p . 



Lemma 12.7( b) and Theorem 13.11 imply that the exact eigenstates of the regularized 
Mattis Hamiltonian are Ur] u n , with the states 77^ in (I2.37p . The corresponding eigenvalues 
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arc 



£n,v = ^ w s(pK(p) + H Q ({q r , s (x)}) + S (3.11) 
»=± P eA s * 

with ifg({g riS (:r)}) as in (13. 9p . but with the operator Q r>s (x) replaced by the eigenvalue 
q r , s {x) = (u +tS (x) + ru- tS (x))/y/2. In particular, the Mattis model has a non-degenerate 
groundstate Ufl with energy ho- 
using these results, one can compute the partition function Zp = Tijr ( e~^ H \ and free 

energy Qp = — \n(Zp)/ (3 of the Mattis model, with (5 the inverse temperature and the 
trace Trj- in J 7 . The additive form of the diagonalized Mattis Hamiltonian in (13 .6p implies 
Zp = Z B ,pZQ t pe~^ £o , with Zx,p = Trx(e~ l3Hx ) for X = B,Q, and with partial boson- and 
zero-mode traces 

Tr B (& B ) = 2^(7? n ,o, B Vn,o) 

~ df A ( 3 - 12 ) 

V 

for boson operators Ob (that only depend on the b^\p)) and zero-mode operators Oq 
(that only depend on the zero-mode operators N r ^ s (x) and R riS (x)), respectively, with r] n „, 
n and u defined in Lemma 12.7( b). Furthermore, this implies that the free energy has 
the form flp = Q Bj p + ^Q,/3 + £o with the groundstate energy So in (I3.10p . One finds by 
straightforward computations (see Appendix IC.3I for details): 

Result 3.2. The boson- and zero-mode contributions to the free energy of the Mattis model 



are 



n 



B,p 



^^lln(l-e-WP)) (3.13) 



Oqj> = -^H^Tj) ~VR ln ^ + °( L3e ~ cL ) ( 3 - 14 ) 

art \Ri,„\/A/ 2B 



for some c > that is independent of L, with w s (p) in ^3. 7|) ; and A and Vf in (EZ3j- 



{Proof: The formula in ( 13. 131) is a simple consequence of a well-known result stated in 
(IC.21I) . The derivation of (13. 14f) can be found in Appendix IC.3.11 ) 

Remark 3.3. Note that the free energy density flp/L 2 has a well-defined IR limit (since the 
factor 1 / L 2 turns the p-sum in (I3.10p and (13 . 13[) to Riemann sums converging to well-defined 
integrals), and the zero-mode contribution VLq^/L 2 vanishes in this limit. In Section HI we 
give simpler formulas for the free energy density obtained in the QFT limit L — > oo and 
a 0+ 



For future reference, we summarize a few useful identities that are simple consequences of 
the above definitions and Theorem I3.1ll"1 To compute thermal expectation values (O) p = 

11 Formulas like the ones given in this paragraph are also used in [20] . 
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Ti>-(e P H 0}/Zp of operators O on the fermion Fock space T one can use the (trivial) 
identity 

(0) p =^-Trjr(e- p& 6), = U^OU (3.15) 

with H in ( I3.6p . Therefore, one only needs to know how to express O in terms of boson- and 
zero-mode operators. Moreover, for operators such that O can be factorized into a boson- 
and zero-mode part Ob and Oq, respectively, the thermal expectation value factorizes 
accordingly: 

P (3.16) 
(Ox)x,p = ^—Tr x (e^ Hx O x ), X = B,Q 

(this is a simple consequence of the definitions, including (I3.12p ). We are particularly 
interested in thermal expectation values of products of "time" evolved density- or fermion 
operators A(t) = e lHt Ae~ lHt with complex t (in applications, t is usually taken to be either 
real or purely imaginary, and our results allow for both). We use the trivial identities 

(Ai(h) ■ ■■A N (t N )}p = l-Trjr(e- pA Mh) • • -A N (t N )), N e N, 

z fi (3.17) 

A(t) = U ] e lHt A e-' lHt U = j tk U ] AU e~ ' lt& , t E C 

implied by (13. 6 p and (13. 15ft . We refer to such thermal expectation values as (thermody- 
namic ) correlation functions, o 

For the computation of correlation functions below, we will need A(t) for the operators 
A = N rtS (x), R rtS (x) and J r , s (p) introduced in Section 12.21 The result is given by the 
following: 

Proposition 3.4. Let r, s = ±, t G C, and U and H as in Theorem \3.1\ Then, for all 
x e Am, 

e iAt U f N r>s (x)Ue- im = N r>s {x) (3.18) 

and 

e 1&t U ] Rr, s (x)U e~ i6i = e irkrAx)t/2 R r , s (x) e irk ^ {x)t/2 e ->rM*t)+«.(*tf] (3.19) 

with 



k r , s (x) = VJ ^Y. [M 1 + Ti) + rr'{l - 1i))N r , tS {x) 



r'=± 



2 

+ 1 E ( T ^-N r , !S (y)+ l2 N r , ! _ s (y)) 

2/6Aid ^ 1 



(3.20) 



and 



12 The extension of our results to include possible time orderings is straightforward. 
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the index set in the last sum is defined as 

A* D ^{peA* D :p^0}. (3.22) 



Moreover, for all p G A*, 

L 



-iw s (p)t 



+ <7 Z PT^(-P)e i ^ (p)t ) - */| T ^-x(p)0+,-.(p.)^_,o 

/ V a 1 + 7i 



(3.23) 



with 



s > ( n • / l rr / J /^f(1 -7iX(p)) . / w fl /(p) , 
V 8vr \ y uv(p) V v F (l -7iX(Pj) 



and 



U s . 



(Proof: See Appendix IC.2I ) 

Note that in the special case 7i = 72 = and t — 0, (I3.23p -f l3.25p reduces to ( 12.26)) . 

3.2 Density correlation functions 

We now consider density correlation functions or, more specifically, the two-point function 

(J ruS1 (x u t 1 )J r2 ,s 2 (x2,t 2 )}p } J r ,s(^t)= e iHt J r ,s(x)e- [Ht (3.26) 

with the densities J rs (x) in (I2.14p . 

Similar to the fermion field operators, the J r , s (x) are operator- valued distributions and 
thus do not automatically have well-defined products. When computing (13.261) . it is there- 
fore convenient to use the regularized densities J r , s (x; e) in (12.141) : for e > 0, the latter are 
operators with well-defined products on the domain T> discussed in the paragraph after 
Lemma O Calculating J rs (x,t;e) = e^W J r , s (x; e)Zie^ using fl2TB"UD . (IZ32]) . ( ET32]) . 
fl3TT5|) . and (JS23D, yields 

J r , a (x, t; e) = J+ (x, t; e) + J~ (x, t; e) + J r ° (x, t; e) 



J (x,t;e) =7^ V V <' s (p)6 s '(p)e 



ipx-kj s /(p)i e -e|p s |/2 



P6A s -'-± (3 _ 27) 

J r ° s (x,t;e)=-^(iV r , s (x_ s )-^ ^ /'(a; a -y;e)JV rVa (y) 

r'=± 2/eAiD 
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with the function 

f{x.e) = ^JL-£ £ e^e-H/ 2 (3.28) 

and superscripts +, — , and indicating the creation-, annihilation-, and zero-mode parts, 
respectively. Note that f'(x; 0+) = 7 2 [4,o - (a/L)]/[2(l + 71)]. 

With that, one can use (I3.17P to compute the correlation function in ( 13 . 26j) . The boson 
part of the thermal expectation values can be computed by a well-known result stated in 
Appendix lC.lt Lemma IC.3( a). The zero- mode contribution is computed in Appendix lC.3.2[ 
We summarize the result as follows. 

Result 3.5. For j = 1,2, let rj,Sj = ±, e A S; and tj G C. Then the density two-point 
function in A3.26\) is given by the limit e — > + of 



1^ \l) H"^(P)C(P) e ^(p)_i +<,si(-P)^ 2 ,s 2 (-P) 
peA^nA*, »=± ' (3^9) 

e -ip(x l -x 2 ) e — e(bail+[Ps 2 l/ 2 _|_ 0(L~ l ) 



g-iw s (p)(ti— 1 2 ) 



I — e -/3oj s (p) 

withu s (p) in fiS/Kj andv^' s (p) in {3.21$. 



As seen in Appendix IC.3.21 the zero-mode contribution to the correlation function 
(J° S1 (xi, t\\ + ) J° 2 S2 (x 2 , t2j + ))q,^ vanishes like 1/L in the IR limit (represented by the 
term 0(L _1 ) in (l3~29l ). Note that the p-sum in (13.291) is a Riemann sum converging to an 
integral in the same limit. 

Since the density operators are linear in the boson operators, it is straightforward to 
compute thermal expectation values of an arbitrary number of "time" evolved density 
operators using Wick's theorem (the zero-mode terms do not contribute in the IR limit; see 
Appendix I C. 3 j) . 



3.3 Fermion correlation functions 

A general Appoint fermion correlation function for the Mattis model is given by 

(^^'^•■•C.ir^^V ^>,t) = e iHt ^)e- iHt (3.30) 

for arbitrary AT e N, with qj,Tj, Sj = ±, x 3 - 6 A s ., and tj 6 C; here and in the following, 
we use the convenient notation [20J 

V£„(x) = <M X )' ^ s (x) = < s (x). (3.31) 

Proposition 12.81 gives regularized fermion field operators ^ s (x, e) = V r)S (x; e)^ 1 in terms of 
the operators J r>s {p) and i2 r>s (x_ s ). This, Proposition 13.41 and straightforward computa- 
tions, yield the following expressions for the fermion fields ^,s( x ; ^ e ) = elHt ^V?,s( x > e )^ e_1 
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q = ±, 



Z $<1 ( x p A e -iq[K% a (x,t;e)+K riS (x.,t;e)] 



^ (X ' tie) = ^ r ' s ^^ (3.32) 



with 



K-(x,t;6) =«.(*_.,*) -ir^ £ £ <M &s , (p)e ipx-^(p)* e -e bs |/2 

K r +(x,t;e)=K-(x,t;e)t (3.33) 
if° s (x,t; e) =r-^\x a N rtS (x- a ) - _ 2/5 e ) N r'-'{y)) ~ k r , s ( x -s)t 

r'=± j/GAid 

and 

f( r - r )- j 72 Q V ip^p-^H/ 2 

/(X,e) - 1 2(l+ 7l )L 2L. p e e (3.34) 

where, again, superscripts +, — , and indicate the creation-, annihilation-, and zero-mode 
parts of an operator K r>s (yL, t; e) (cf. (I3.27P ). Note that the operator J r s (x, t;e) in (I3.27p . 
multiplied by a factor 2irar, is equal to the partial a^-derivative of -fT riS (x, t; e). 

With that, (13.1 7p provides a convenient recipe to compute the fermion Green's function 
in (I3.30p : by (I3.16P and (I3.32p . one has to compute the boson expectation value 

and the zero-mode expectation value 

(S* tl (x 1; tr, e) • • • S% tSIf (x N , t N ; e)) Q)fj . (3.36) 

Since K^ s (-x,t;e) is a linear combination of boson operators, the boson part in (I3.35P 
can be computed using well-known results stated in Appendix IC.14 Lemma IC.3( b). We 
note that all terms involving u s (x_ s ,t) are 0(L _1 ) and thus do not contribute in the IR 
limit. 

The zero-mode part in (I3.36P can be obtained by straightforward computations; see 
Appendix IC.3.21 for details. One finds that, in the IR limit, this term simplifies to 

(o, R rusi (x^ si y^ ■ ■ ■ R rN)SN (x NjS _ N y» r »n). (3.37) 

We thus obtain the following: 

Result 3.6. For N e N and j = 1, 2, . . . , N, let qj, rj, Sj = ±, Xj G A s , and tj 6 C. Then 
the fermion correlation function in h3. 30\) is given by the limit e — > + of 



3=1 

x ( II G ri,s M (xj ~ x fe , tj - t k ; e)"**) (l + 0(L- X )) 



(3.38) 



l<j<k<N 

22 



where 

In a i)Slj r 2 , S2 (x, t- e) =nr 2 a ^ (j) E^rk, S i(P)C(P) 



c 



iw s (p)t p -iw a (p)i 



"iP x p - e (|Psil + |Ps 2 l)/ 2 



P>s (e) = 27rSea, s ,r,s(O,0;e) (3.40) 
wift a; s (p) m pT^ ) and v*' s (p) m flff.^l ). 

Note that the p-sum in f !3.39|) is a Riemann sum converging to an integral in the IR limit. 
Even though this integral is singular, it is well-defined since its definition also specifies how 
to treat the si ng ularity; see Section ID. 2. II The terms 0(L _1 ) in ( I3.38|) are mainly zero-mode 
contributions^ and they are computed in Appendix 10.31 

We emphasize that, in the IR limit, the effect of the zero-mode terms is negligible except 
for the factor in (I3.37p . which is independent of temperature and coupling. This factor can 
be easily computed by "moving Klein factors to the right" using 




XrA*)Xr>A* ) ~ <! -rV , ^ ™ otherwise [3A ) 



(this follows from d2.34[) ) and (12 .43 p . Thus, a fermion Green's function in (I3.38j) is non-zero 
only if N is even and if, by using (13.4 ip repeatedly, the product of Klein factors in (I3.37P 
can be simplified to either 1 or —1. In particular] 14 ! 

(fi, R riSl (Xx) ' -R r2 S2 (^2) ' fi) fiqi,— q2^r\,r2^si,S2^xi,X2 (3.42) 

<n,i2 ri , ai (a; 1 ) ri « 1 ...i2 r4ja4 ( a ;4) r454 n) = 

"91,— g2"ri,r2"«i, s 2 "a;i, x 2 ^qs, -q 4 ^r 3 ,r 4 ^s 3 , s 4 $x 3 ,x 4 ^ 

^<?i i — 93 "»1 ,^3 "ai ,S3 &xi ,X3 3q2,-q 4 <V 2 ,r 4 <^s 2 ,S4 "x 2 ,x 4 
"I - 0q\ ,—q 4 $r\ ,r 4 &si ,s 4 ^xi,x 4 ^q2,—q3^r2,r^^S2 ,S3^xs,X3 ■ 



4 Quantum field theory limit 

The results of the previous section allow us to compute, in principle, any quantity of physical 
interest for the Mattis model. However, the given formulas are in general complicated, and 
it is not easy to discern from them the qualitative features of the model. For this reason, it 
is useful to only consider the physical behavior for length scales that are much larger than 
the UV cutoff a, but still much smaller than the IR cutoff L. Indeed, as will be shown in 
this section, one can obtain simplified formulas in the QFT limit L — > oo and a — > + (after 
suitable multiplicative renormalizations). As illustrative examples, we analyze in detail the 
free energy and fermion two-point function of the Mattis model. Computational details can 
be found in Appendix IDl 

13 The other contribution comes from u s (x- s ,t), as already mentioned. 
14 Below we write Xj short for Xj t - Sj . 
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4.1 Free energy 



Computing the contributions to the free energy in leading order in a and l/L we obtain: 
Result 4.1. The QFT limit of the (renormalized) free energy contribution ttp — So = 

Q B ,P + &Q,p is 

lim lim hL~ 2 (to* - S ) = (4.1) 

a^O+i-Kx, 3v F VAfi 2 

with the parameters A and v F in h3.8\) . 

{Proof: See Appendix ID. II ) 

This result is easy to derive when 72 = (A = 1) and, in this case, the corrections to 
the leading term of aL~ 2 (Qp — So) (at non-zero a) are exponentially suppressed (i.e. they 
are O(e-^ /a ))0 

As seen in the proof in Appendix ID. II the result for non-zero 72 can be obtained by 
replacing the original dispersion relation (13. 7\\ by the effective dispersion 

w±(p) = v F VA\p±\ ( 72 x(p)p+P- ^ 0) (4.2) 



(corresponding to renormalizing 71 — > yl — A (1 — 7^) and 72 — > in (13. 7p ) in the expres- 
sion ( 13.131) for the boson free energy. This can be understood as follows. The dominating 
contributions to (13.131) come from terms for which the dispersion is minimized. A glance 
at Figure [2] shows that, in this regime, the dispersion relation is excellently approximated 
by 63). 



4.2 Fermion two-point functions 

As discussed above, we are interested in the behavior of correlation functions for interme- 
diate length and time scales defined by the regime 

a < ||x| ± v F t\ < L, a < ||x| ± v F y/At\ < L. (4.3) 

In this section, we will focus solely on the case of fermion two-point functions in this regime 
and at zero temperature. 

To illustrate important features of the QFT limit, we first consider the special case 
72 = for which simple formulas can be obtained even before taking the limit a — > + . We 
will need the following special function 



a(z) = e~ E ^ z \ 



/oo 
(l/t)e- zt dt (arg(^)<7r), (4.4) 



15 As seen in [3], the Fermi velocity vf is proportional to a when the Mattis model is derived from the 
original lattice fermion model. Thus, in this case, the corrections to (|4.ip do not vanish in the UV limit 
a — > + , and Result |4~T1 should be regarded as the low-temperature limit of the free energy. 
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whose asymptotic behavior is given by (see Appendix ID. 2. lj) 



a{z) = e^z(l + 0(z)), a(z) = 1 - — (l + O^ 1 )), (4.5) 

z 

with Euler's constant 7 = 0.5772(1). 

Result 4.2. Let 72 = ; qj, rj, Sj = ± for j = 1, 2, x G IR 2 , and t G C. Then the fermion 
two-point functions at zero temperature and in the IR limit are 

lim (^l, Sl (^,t)ipfl S2 (0,0)) = S qi - q2 6 rur2 5 Sl , S2 ^ x _ sii0 (e 7 7r) 1_A 



. . 1 1 / a 



,2 \ (K-l)/2 (4.6) 



2tt 0+ - i (nx Sl - v F t) V(0+ + i£ F t) 2 + (x Sl ) 2 
with vf in /IS. 8\) . 

*=i(* + iT>), B = ^f, (4-7) 

and 

F , r ,x _ "(f (° + ~ - M))) (y+1)/ V(f (0 + + i(rx + n^))) ( ^ 1)/2 

-(f(0+-i(ra;-^))) " ^ 



(Proof: See Appendix ID. 2. II ) 

It is instructive to compare the short- versus the intermediate length- and time scale 
structure of the correlation functions in ( 14. 6p . Using the first relation in ( 14. 5 p one finds for 
short length- and time scales 

lim <^(x,t)^ s (0,0)> oo =^_ s ,o; 1 ' 



L^oo^ r > sX ' JTr > sX ' "™ a x - s ' u 2vr0+ -i(rx s - v F t) (4.9) 

x (1 + 0(x s /a) + 0(v F t/a) + 0(v F t/a)) , 

identical with the short- distance behavior in the non-interacting case. 

From the second relation in ( 14. 5ft . it follows that F r>s (x,t) — > 1 in the regime H4.3f) and 
in the QFT limit. The result in ( I4.6P thus shows that the fermion two-point functions in 
the Mattis model, for 72 = and in the QFT limit, have an algebraic decay with exponents 
depending on the coupling parameter 71. We note that such behavior is one of the hallmarks 
of Luttinger-liquid behavior [9]. 

Note that K > 1 for the interacting case, and therefore, to get a non-trivial UV limit 
in the regime (I4.3p . one needs to multiply (14. 6 j) by a factor oc a^~ K ^ before taking a — > + . 
We find it convenient to choose this multiplicative renormalization factor as (e 1 7rLo/a) K ~ 1 
with a finite, arbitrary length scale parameter Lq. 

Remark 4.3. The computation above provides a generalization of the limit in which the ID 
Luttinger model reduces to the massless Thirring model; see e.g. 



25 



One of our main results is that, after a multiplicative renormalization similar to the one 
above, such an algebraic decay of fermion two-point functions holds true in the QFT limit 
even for non-zero 72: 

Result 4.4. Let qj,Tj,Sj = ± for j — 1,2, x 6 M. 2 , and t e C. Then the (renormalized) 
fermion two-point functions at zero temperature and in the QFT limit are 

Jim lim — - — (x, f)^ S2 (0, 0)> = 5 9l _ 92 5 rijr2 5 SljS2 6(:r_ sl ) 

1 1 / ^ \^-D/ 2 (4.10) 

xC(7i, 72 )- 1 



2tt 0+ - i(nx Sl - v F VAt) \(0+ + kVv / It) 2 + (x Sl ) 2 / 
with the dimensionless constant (7(71,72) given in Appendix \D.2.S\ 

K^J^+S), B = J*E*. (4.1i; 
2VS Vl +7i 



and A, m A3. 8\) : L > is an arbitrary length scale parameter. 
(Proof: See Appendix ) 



The constant C(7i,7 2 ) is defined by an integral of trigonometric functions (see (ID. 251) - 
(1D.26j) ). The special case 7 = 71 = 72 is plotted in Figure [3j 

The similarity between the correlation function in ( 14.101) and the corresponding result 
one would obtain for 72 = (A = 1), together with the result derived for the free energy 
in section |4~T1 might leave the reader with the impression that, in the QFT limit, no new 
physics is obtained by allowing 72 7^ 0. We stress that this is not so. It is true that, 
apart from the overall multiplicative factor (7(71,72), Result fl~4l can be obtained from the 
corresponding result for 72 = by effectively replacing v F by vp\fA and K in ( 14.7)) by that 
in (14. lip . However, it is not possible to find an effective coupling parameter 71 such that 
the result in f 14 .10)) is obtained from the one for 72 = and by changing 71 to 71: the UV 
limit of the Mattis model depends non-trivially on both coupling parameters 71 2- 



4.3 Other correlation functions 

Our results in the previous section suggest that the following renormalized fermion operators 

*t(x) = (e^Lo/a^K-W^ix), (4.12) 

with K in (14. lip and L > an arbitrary length scale parameter, have well-defined non- 
trivial correlation functions in the QFT limit: 

Conjecture 4.5. Let qj, Vj, Sj = ±, x, G M 2 ; and tj G C for j = 1, 2, . . . , N, and N G N. 
Then the following renormalized fermion correlation functions of the Mattis model in the 
QFT limit, 

lim lim (^(xLtj) ■ • .*«^(x*,ttf)> (4.13) 
are well-defined distributions for all 1//3 > 0. 
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0.75 




Figure 3: The dimensionless constant (7(7,7), occurring in the QFT limit of the fermion 
two-point functions, for the special case of 7 = 71 = 72- 

A proof of this conjecture, which would imply that the (constructive) QFT limit of the 
Mattis model exists and is non-trivial, involves taking the limit a — > + in Result 13.61 for 
arbitrary N. The calculation might be done in a similar way as the N = 2 case presented in 
the previous section and Appendix ID. 21 However, due in part to the increasing complexity 
in the contribution from the Klein factors to higher order correlation functions (see for 
example Equation (I3.43[) corresponding to iV = 4), a formula for arbitrary iV will not be 
attempted in the present paper. 

It is also possible to compute the QFT limit of density correlation functions. Here we 
only give the result for the case 72 = and in the IR limit (see also the remark at the end 
of Section 13. 2p : 

Result 4.6. Let 72 = 0, 77, Sj = ± forj = 1, 2, x e IR 2 ; and t £ C. Then the (renormalized) 
density two-point functions at zero temperature and in the IR limit are 



2 



Km a<J ri)51 (x, t)^a(0»0)>oo = S n,*^ S *-'ifl777ZA 

I B + r 1 r 2 B~ 1 + (77 + r 2 ) B + r 1 r 2 B~ 1 - (77 + r 2 ) 
X I (0+-i(x Sl -v F t)) 2 + (0+ + i(x Sl +v F t)) 2 



(4.14) 



with 



e r y(x,t) = (^r^j (^A5 r yai (0 + — i (rx — v F t)) 

((B + rr'B" 1 ) + f (r + r')) a x Q (0 + - i {fx - v F t)) 



E 



(4.15) 



the function ol\{z) = z 2 e z (l + z), and the parameters v F and B given in ( Iff. 8\) and 
( l^.iip respectively. 
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(Proof: See Appendix ID. 31 ) 



The corresponding QFT limit of (14.14p should follow in a direct manner (we expect that 
e r r /(x,t) — > in a distributional sense as a — > + ). Finding a simple analytical expression 
of the density two-point functions for general 72, and in the full QFT limit, appears more 
complicated and will not be attempted here. 

4.4 Positivity 

As mentioned, the massless Thirring model can be obtained from the ID Luttinger model 
as a limit that is similar to the full QFT limit of the Mattis model. One can therefore regard 
the model obtained from the Mattis model in the full QFT limit (referred to as Mattis QFT 
for short in the following) as a 2+1 dimensional variant of the massless Thirring model. 
The Thirring model has been studied extensively in the context of axiomatic QFT (see e.g. 
[34J and references therein), and it is therefore interesting to also discuss Mattis QFT from 
this point of view. 

Obviously, Mattis QFT is not Lorentz invariant, and it is therefore not a QFT theory in 
the sense of axiomatic QFT@ For this reason, many issues of axiomatic QFT do not apply 
to Mattis QFT. One important issue that does apply, however, is positivity (in the sense of 
Osterwalder and Schrader [35J ) . We now sketch why positivity holds true for Mattis QFT. 

Our results in the present paper show that the Mattis model is defined by a self-adjoint 
Hamiltonian on a separable Hilbert space, and this Hamiltonian is bounded from below. 
Moreover, the fermion field operators are well-defined operator valued distributions on this 
Hilbert space, i.e., smearing them by suitable testfunctions gives bounded operators^ For 
these reasons, positivity is obvious for the Mattis model. Since positivity is preserved under 
limits (see e.g. [34J), the result follows. 

Our argument above is sketchy since the details are very similar to the ones in the proofs 
of positivity of the massless Thirring model in [331 134"] . 

5 Final remarks 

1. In this paper, we presented a mathematically precise construction of the Mattis model 
formally defined by the regularized Hamiltonian in (jl.2p . As proposed in [31 0], the Mattis 
model provides an effective description of 2D lattice fermions in a partially gapped phase 
close to half filling. We diagonalized this Hamiltonian, and we computed all thermodynamic 
correlation functions. We also proved that the fermion two-point functions have a non- 
trivial QFT limit. 

16 Note that axiomatic QFT has been mainly concerned with Lorentz invariant QFT models motivated 
by particle physics — the fact that the full QFT limit of the ID Luttinger model yields a Lorentz invariant 
QFT model is a coincidence. 

17 We do not use test functions to deal with the distributional nature of the field operators but another 
method similar to the one in |34j . 
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One of our aims was to develop computational tools that can be straightforwardly 
generalized to more complicated models of similar type. One such model is a spinfull 
variant of the Mattis model that can be derived from the 2D Hubbard model by a partial 
continuum limit Fj 



2. Although the Mattis model is similar to models of parallel Luttinger chains, there is one 
important difference: the Mattis model describes two sets of parallel Luttinger chains that 
are orthogonal to each other. Apart from the limiting case 72 = (for which the Mattis 
model is essentially equivalent to the model studied in [27]), the two sets are non-trivially 
coupled by density- density interactions. 

3. As mentioned in the previous remark, the orthogonal Luttinger chains of the Mattis 
model are not coupled for 72 = 0. This explains the locality of the fermion correlation 
functions parallel with the Fermi surface directions in this case. We find that this locality 
holds true even at non-zero values of 72. This is due to the 8 Xj jXk -factors that appear in 
the Klein factor contributions in (13.421) and (I3.43p . For density correlation functions, this 
locality seems to be restricted to the case 72 = 0, at least for length scales of order a. We 
suspect that this locality of the fermion correlations functions can be changed by adding 
correction terms to the Mattis Hamiltonian (see below). 

4. The results obtained in the present paper allow, in principle, to check the validity of 
the approximations used to obtain the 2D Luttinger model from the 2D t-t'-V model [3] 
(similarly as done in [15] for the model proposed by Luther [H]) and, if necessary, to improve 
upon these approximations. To be more specific we discuss one important approximation 
used in [3], namely replacing the full band relation by linearized ones (the approximation 
in (IA.8j) ). As elaborated in Appendix |El the leading corrections to this approximation are 
the following terms that can be added to the Mattis Hamiltonian, 



r,s=± Js 



(5.1) 



and 



H 3 = - ^ v F d 2 x : (a 0)2 aV4 lS (d-A,J ~ a i,2a 2 ripl s (id s )(d 2 s ip rs ) J : 

r,-=± Js (5.2) 

(<9?. s Vv |S )(x) = (Vr,s(x + ae_ s ) + ip r>s (x - ae_ s ) - 2ip r>s (x) j 



with dimensionless constants atjk given in flE.2j) . The correction terms in ( 15. ip and ( 15. 2 p 



take into account curvature effects of the band relations orthogonal to, and parallel with, 
the Fermi surface, respectively, up to terms 0(a 4 k 4 ). The former can be bosonized as 
follows (the interested reader can find details in Appendix IE.2p . 



+ 



a 3>0 a 2 



i (5.3) 



x (n 4 + 6H^ S $J 2 + (<9 S $J 4 ) + (<9 S 2 $ S ) 2 + (dA) 2 ] + 0(L 



-1\ \ x 



18 J. de Woul and E. Langmann (in preparation). 
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(the corrections are zero- mode terms and terms given in Appendix [E]). Thus H2 corresponds 
to non-linear boson interactions. Note that the coupling constants in ( 15. 3p are of higher 
order in the UV cutoff a than in (15. f p . This suggests that the corrections H2 are less 
important than one might naively expect. 

It is also possible to bosonize H3 in (I5.2p using Proposition 12. 8[ but for this term the 
Klein factors R r)S (x- s ) do not cancel. We believe that these Klein factor terms are important 
corrections to study since they can make fermion correlation functions non-local parallel 
with the Fermi surface directions. 

In a similar way, one can bosonize correction terms for all other approximations that 
lead from the t-t'-V -model to the 2D Luttinger model (see Appendix|A]for more details on 
these approximations). 

5. The term proportional to ao,2 i n H 3 describes hopping between nearest neighbor chains. 
Previous work on coupled Luttinger chains suggests that Luttinger-liquid behavior can be 
destroyed by arbitrarily small interchain hopping; see e.g. [37J[38j[39] and references therein. 
However, in our case the hopping parameter vanishes in the UV limit a — > + , and it is 
therefore an open question if H 3 can change the qualitative behavior of the Mattis model 
in the UV limit. 

Note that there is a particular location of the Fermi surface of the 2D t-t'-V model where 
a o,2 — (see Appendix IE. ip . In this case the leading correction in H 3 to the Mattis model 
is zero. 

6. The Hamiltonian methods used in this paper have several merits, but renormalization 
group methods based on functional integrals [lOj HI] seem more powerful when it comes 
to determining the effect of correction terms that spoil exact solubility (like the ones in 
(15.ip - (j5.2p ). We believe that it would be fruitful to try to combine these methods with 
existing functional integral approaches to bosonization [121 H3j HI] . 

7. Anderson in an influential paper [TO] (appearing three years after Mattis' [lj) suggested 
that 2D fermions can have Luttinger-liquid behavior. Our Result |4~41 shows that the fermion 
two-point functions in the QFT limit have a non-trivial algebraic decay. This property is 
often regarded as a hallmark of a ID Luttinger liquid [9]. However, 2D Luttinger-liquid 
behavior is subtle (one difficulty being that different definitions seem to be used in the 
literature) . 

A mathematical criterion for Luttinger-liquid behavior in 2D in the context of construc- 
tive QFT was given by Salmhofer [45j. Since our methods are different from the ones used 
in constructive QFT, it is not easy to directly apply the criterion in |45j to our results. 
However, Mastropietro [27] has shown Luttinger-liquid behavior in the sense of [H] in a 
model which is, essentially, the limiting case 72 = of the Mattis model. His result, to- 
gether with our Result I4.4[ suggest that the Mattis model has such behavior for all allowed 
values of the coupling parameters 7^2. This is consistent with earlier results suggesting 
that non- Fermi-liquid behavior is only possible in 2D models with Fermi surfaces that have 
flat parts; see [4HI |4"7] and references therein (note that the latter results are restricted to 
weakly coupled fermion systems). 
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We emphasize that, even if one proves Luttinger-liquid behavior for the Mattis model, 
this does not prove Luttinger-liquid behavior of the 2D t-t'-V model. For that one would 
need to substantiate the physical arguments in [31 H] by mathematical proofs. 

8. The most general form of the Mattis Hamiltonian that can be solved exactly using the 
methods developed in the present paper is 

H = v ? E E(x) 2r ^ : ^ (k) ^ (k): 

r,s=± keA; 

! 2 . (5.4) 

+ Y (jj ^WV(P) : ^r, a (-p)^r', S '(p): 

with suitable interaction potentials V r s r i ^(p) = K-',s',r,s( — p) (the notation here is the same 
as in Definition 12. ip . It is straightforward to write H in (j5.4p as a non- interacting boson 
Hamiltonian, and by diagonalizing the latter one finds necessary and sufficient conditions 
on Vr,s,r',s'(p) such that the generalized Mattis model is well-defined (these conditions are 
similar to those for the Luttinger model; see e.g. |20j, Equations (3.22)-(3.24)). 

The model discussed by Mattis [1] and Hlubina p3] correspond to the special case 
K-,s,r',s'(p) = ^(p) (independent of r,s,r',s'), with the latter author also assuming that 
the interaction is long-ranged, i.e. V(p) = for |p| > q Q and go *C n/a. To not overburden 
the presentation, and due to the motivation provided by our previous work [31 H], we 
restricted our discussion in this paper to the case for which the interaction is short-ranged, 
i.e. K-.s.r'.s'Cp) * s essentially independent of p. It can be expected that the Mattis model 
with long-range interactions has physical behavior that is qualitatively different. 
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thank Farrokh Atai for carefully reading the manuscript. 



A Relation to lattice fermions 

As discussed in the introduction, one motivation for the Mattis model is that it can provide 
an effective description of a 2D lattice fermion system, as proposed in [21 H]. Here we 
describe in more detail the lattice fermion system and the nature of the arguments on 
which our proposal is based. 

We first introduce some useful terminology: we define fermion models by giving the 
so-called set of one-particle quantum numbers S, the normalization of the CAR (canon- 
ical anticommutator relations) z > 0, the dispersion relation and interaction vertex 
VkiMMM- By this we mean a fermion model defined by fermion field operators ip k , k G S, 
obeying CAR with {V'fcjVv} = z5k,k>, and with the Hamiltonian 

kes kiMMM&S 
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We note that this defines the model completely if (and only if) S is a finite set. 

The lattice model describes spinless fermions on a square lattice with nearest-neighbor 
(nn) and next-nearest neighbor hopping and nn density- density interactions. This so-called 
t-t'-V -model can be defined as follows: The set of one-particle quantum numbers is equal 
to the Brillouin zone 



BZ d = f <jk=(fc 1? A; 2 ): k ± e j- (z + ~Y ~ < * lf2 < ^ (A.2) 



of the lattice, with k± = (k\ ± A^V^; here a > is the lattice constant and L > the 
lattice size such that L/(2a/2o) is a positive integer. Choosing the normalization of the 
CAR as \Lj (27r)] 2 , the band relation of the model is 

e(k) = — 2t[cos(a/ci) + cos(afc 2 )] — 4t' cos(a/ci) cos(afc 2 ), (A. 3) 

and the interaction vertex is 

«(ki, k 2 , k 3 , k 4 ) = M(ki - k 2 ) ^ \Y) ^ki-^+ks-k^Trn/a (A.4) 

with -u(p) = a 2 V[cos(api) + cos(ap 2 )]/(87r 2 ) the suitably scaled Fourier transform of the nn 
interaction potential; see Section 4.1 in [3] for further details and explanations. The box 
surrounding Figure Q] shows the boundary of the Brillouin zone in (IA.2I) . with the horizontal 
axis corresponding to k\ and the vertical axis to /c 2 . The vectors k are interpreted as 
momenta of the fermions. 

One important parameter of the t-t'-V -model is filling u, which is defined as the average 
number of fermions per lattice site; it is in the range < v < 1. Various arguments 
(discussed in [3]) suggest that, at half-filling (i.e. for v = 1/2) and for sufficiently large 
values of V/t, this model describes an insulator, i.e., all fermion degrees of freedom are 
gapped. The challenge is to understand the model away from, but close to, half-filling. 
Experimental and theoretical results (discussed in [3]) suggest that fermions associated with 
different regions in the Brillouin zone can have different physical properties for v ^ 1/2 
and, in particular, fermions close to points 

Qrs = (rQ/a,rsQ/a) } r,s = ± (A.5) 

with Q m 7r/2 (the four points that are midpoints of the dashed line segments in Figured]) 
are gapless, while fermions in other regions can still be gapped. This suggests to rewrite the 
fermion model by dividing the Brillouin zone in eight different subsets Q r s + A* s , r = ± 
and s = 0, ±, 2, with Q r s conveniently chosen points such that 

E/( k ) = E E /(Q^+ k ) ( A -6) 

keBZ r,s keA* s 

for any function / on the Brillouin zone; see Figure 2 in [3] for how this division is defined. 



In particular, the regions A* s , r,s — ±, associated with the points in (1A.5I) are rectangles 
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defined by —-nja < k^ s < ir/a, and a similar constraint in the /c s -direction (see [3], Equa- 
tion (42)) with a = \/2a/(l — k) and < k < 1. This parameter k determines the size 
of the regions treated by bosonization and it also gives the UV cutoff of the Mattis model 
that we use in the present paper 0, H] ■ 

This division of the Brillouin zone allows to rewrite the t-t'-V -model model as a model 
of eight flavors of fermions, distinguished by labels r = ± and s = 0, ±, 2, with each of them 
labeled by momenta k in different sets A* s : the dispersion relations and the interaction 

vertices of these fermions are e r>s (k) = e(Q r s + k) and 



(ki, . . . , k 4 ) = u(Q ri , Sl - Qr 2 , S2 + ki - k 2 ) 

X ^Qr 1 ,s 1 -Qr2,S2+Qr3,s 3 -Qr4,34+kl-k2+k3-k4,27r/an- 



(A.7) 



We refer to the fermions with r, s = ± as nodal. The fermions with r = ±, s = correspond 
to regions around the points Q +i o = (7r/a, 0) and Q_o = (0,7r/a) and are referred to as 
antinodal; the fermions with r = ±, s = 2 are included for completeness, but we expect that 
they can be safely ignored in the parameter regime of interest to us [3j. The parameter Q 
is chosen such that the points Q riS lie on the Fermi arcs, and one expects that, for v = 1/2, 

q=tt/2 mm. 

To obtain a model that is better amenable to computations we proposed a series of 
approximations [3]: First, to Taylor expand the dispersion relations e r>s (k), r, s = ±, and 
keep only the leading non-trivial terms, e.g., to replace 

e r;S (k) = e(Q r , s + k) -» e(Q r>s ) + rv F k s , r,s = ± (A.8) 

with vf in [3], Equation (6) (Approximation Al in [3]). Second, to replace the interaction 
vertices of the model by 



( L \ 



2 



V 



ri,si,...,r4,S4 



V2W (A.9) 

nez 2 

(A2 in [5]). Third, to ignore the fermion degrees of freedom with r = ±, s = 2 (A3' in [3]). 
One key result in [3] is that, if Q ^ 7r/2, the approximations above lead to a model with 
density-density interactions only, i.e., all interactions can be written in terms of fermion 
densities 

PvAp) = ^L(ki)^ r , s (k 2 )5 kl+P5k2 . (A.10) 

ki,k 2 eA* ^ ' 



r.s 



To obtain a model where the nodal fermions can be treated in a simple way by bosonization 
two further approximations are needed: First, to drop the restriction on k s in the nodal 
sets, i.e., to replace in (lA.lOj) and in the free part of the Hamiltonian 



a*^a: (a.h; 
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with the latter sets defined in ( 12. ip . which can be done after normal ordering (A3 and A4 
in [3]). Second, to replace the (normal ordered) nodal densities in ( lA.lOp thus obtained by 
the ones in (12. 6p : 

:p r , s (p):^ J r , s (p), r,s = ±, (A.12) 

i.e., to insert integer sums (A5 in [3]). The former of these approximations leads to a QFT, 
i.e., it replaces a model with a finite number of degrees of freedom by one with infinitely 
many ones. The model thus obtained was called 2D Luttinger model in [3] since, first, it 
is related to 2D lattice fermions in a similar way as the ID Luttinger model is related to 
a corresponding ID lattice fermion system and, second, in this model the nodal fermion 
degrees of freedom can be treated using bosonization methods. In pE] we proposed to treat 
the antinodal fermions in the 2D Luttinger model by mean field theory, and we presented 
mean field results showing that there is a significant parameter regime where the antinodal 
fermions are gapped. Moreover, we showed that the latter result is robust, i.e., insensitive to 
changes in k and Q. We also proposed that, in this partially gapped phase, the low energy 
properties of the 2D Luttinger model is described (essentially 19 !) by the Mattis model. 

We note that all approximations above either only affect fermion degrees of freedom 
that are far away (in energy) from an assumed underlying Fermi surface (containing the 
points in OA.Sp ). or they amount to changing the Hamiltonian by terms that are formally 



sub-leading in the lattice constant a. Physical arguments suggest that such changes do 
not (much) affect the low energy properties of the model. However, the validity of these 
arguments is far from obvious. It therefore is important to note that the results in the 
present paper provide a starting point to investigate the validity of the approximations in 
greater depth. In Section Remark 4 we explain this for the approximation in ( ]A.8j) . This 
is only one example, and the other approximations above certainly also deserve attention. 
In particular, one might be concerned about the approximation in (|A.9[) : this amounts to 



ignoring, in particular, the fermion interaction terms associated with the interaction vertex 



X 



(u{Q ri , Sl - Qr 2 ,s 2 + k i - k 2 ) - u(Q ri , Sl - Q 

(IlW sr, (A - 13) 

^27T J kl ~ k2 + k 3-k 4 ,0 / J °Q ri , Sl -Q r2 , 32 +Q r3 , sa -Q r4!S4 ,27rn/a 



for Tj, Sj = ±. While these terms are formally sub- leading in a, there are certainly scattering 
processes with momenta ki,k2 close to the underlying Fermi surface but where |ki — ki| 
is of order ir/a. Furthermore, the approximation in ( 1A.12I) amounts to adding umklapp 
processes that are not present in the lattice model [3], and the role of these should be 
further investigated. 



B Boson- fermion correspondence 

In this appendix we summarize well-known mathematical results on constructive bosoniza- 
tion in one dimension, and we show how to obtain from these the results in Section 12.21 

19 With 71 = 72 and more general cutoff functions; see Remark 12.41 
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B.l Bosonization in ID 



In this appendix, we write r, r' — ± for chirality indices, A, A' G X for flavor indices with 
I some index set to be specified later, and k, k! G (27r/L)(Z + 1/2) for ID Fourier modes. 
We consider fermion operators c^ A (k) defined on a fermion Fock space J 7 with normalized 
vacuum state Q, satisfying the usual relations of ID relativistic fermions: 

{c rA {k),cl, A ,(k')} = 5 r y6 A ,A'Sk,k', {c rA {k),c r , A ,(k')} = (B.l) 

and 

c rA (k)tt = cl >A (-k)n = Wk such that rk > 0. (B.2) 
We also introduce so-called densities 

for p G (27r/L)Z; here, and in the following, the colons indicate normal ordering as in (\2A\i . 

Proposition B.l. (a) The densities in W.3\) are well-defined operators on T and satisfy 
the commutator relations 

Lp 

L?Vl(p), jr',A'(p')] = rS r y8 AtA '— $p+p>,0- 

Moreover, j r ,A(pY — Jt,a{—p) an d 

j rj A{p)fl = Wp such that rp > 0. 

(b) The following operator identity holds true 

rk :cl A {k)c r A {k): = j ^ ■3rA-p)jrAp) ■ 
where both sides are well-defined operators on J 7 . 

(c) There exist unitary operators R t .a on J 7 satisfying the following relations, 

\jr,A{p)i Rr',A'] — r O~r,r'0~A,A'0~p,oRr,A, {R r ,A> ^r',A' } = ^r,r'^A,A' (B-4) 

and 

(n, n n <a a v) = n n <w. ° ( b - 5 ) 

r=± AeX ' r=± A&X 

for all m TiA G Z such that Y2 r A \ m r,A\ < oo. 

(d) The operators j r ,A(~ r P) an< ^ Rr,A, with p G (27r/L)N, r = ±, and A el, generate the 
full fermion Fock space T from fl, i.e. the states 

nnn ] f n n <A n ^ 

K r=± AeX pg^N / V=±AeX / 
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with m rj A £ Z and n r ^A(p) £ No such that J2 r A m r,A < 00 an< ^ J2 P \p\ n r,A(p) < oo ; provide 
a complete orthogonal basis in the fermion Fock space J 7 . 

(e) For e > and -L/2 <x< L/2, 

V rA {x] e) = ^Le r ™ i ^ (0)/i (R r ,A)~ r e rnix ^ A{0)/L 
V27re 

xexp(r^ £ M^ e ^ ^ 
P e^z\{o} P 

define unitary operators on J 7 . Moreover, 

1 /" L/2 



c riA (fc) = lim — / dxV r>A (x; e) e- ifcx (B.8 



L/2 



where the limit is understood in the strong operator topology (see e.g. l3Uf). 



Proof: Proofs of (a)-(d) based on vertex operator algebras can be found in [25], Sections 5.1 
and 5.2, for example (the details of the proof given there are only in the special case of a 
single fermion species, but the generalization to our case is trivial). An algebraic version of 
(e) is also stated and proved in [25], Sections 5.2, but we find it convenient to use a different 
formulation proved in |23j, for example. □ 

Even though the results summarized in Proposition IB . 1 1 are well-known, the mathemat- 
ical literature in which rigorous proofs are provided is, to our opinion, not easily accessible 
for non-experts. We therefore give below a few remarks, where the first four and the last 
concern the mathematics- and the physics literature, respectively. 

Remark B.2. Parts of Proposition IB . 1 1 are special cases of mathematical results stated and 
proved in [221 121], for example. To be specific: (a) can be obtained from a wedge represen- 
tation of an affine Kac-Moody algebra, (b) is a special case of the Sugawara construction, 
and (c)-(e) can be deduced from standard vertex operator constructions. 

Remark B.3. Similar to what was discussed in Remark I2.2| issues related to the unbound- 
edness of the operators appearing in Proposition IB . 1 1 can be dealt with using general math- 
ematical results proved in [32J ; see also Section 12.31 However, we will not elaborate on these 
results here since overly emphasizing this unboundedness issue would be misleading: there 
exists an alternative rigorous approach based on vertex operator algebras that is purely 
algebraic [25] . 

Remark B.4. A delicate mathematical issue is the distributional nature of the fermion field 
operators, i.e., Part (e) of Proposition lB.il As discussed, our approach to this problem fol- 
lows [23]: we consider regularized fermion operators, proportional to unitary operators, and 
which depend on a parameter e > such that, when e — » + , they converge to the fermion 
field operators in a known distributional sense [23] . We find it convenient to use these reg- 
ularized (unitary) operators when computing correlation functions, i.e. inside expectation 
values, since then the limit e — > + can be done in a fully controlled manner. This allows us 
to suppress some analytical details and thus simplify our presentation. Furthermore, overly 
emphasizing these details would be misleading since the essence of (e) is algebraic rather 
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than analytic: as mentioned in Remark IB.3| there exists an alternative rigorous approach 
based on vertex operator algebras [25J. 

Remark B.5. It is well-known that the Hilbert space obtained using just the boson operators 
is a proper subspace of the fermion Fock space, and to get equivalence one has to include 
different sectors in which the fermions are filled to different levels; see e.g. [IS], end of 
Section I, for a discussion of this point. Since fermion operators do not leave sectors 
invariant, while boson operators do, it is impossible to represent fermion operators in terms 
of bosons only. The Klein factors provide the remedy to this problem; see [25] , Equations 
(5.2.1)-(5.2.4) where the Klein factor is denoted as u. As discussed in the introduction, 
in the physics literature, Klein factors can sometimes be ignored (and often are) in ID 
problems, but this is not true in higher dimensions. A constructive approach to Klein 
factors is discussed in [21], Section 3.4, for example. 

Remark B.6. A standard reference in the condensed matter literature in which the results 
above are derived is [21] • For the convenience of the reader, we therefore give the dictionary 
between our notation and that used in [21] ("a -H- 6" below means "the symbol a in [21] 
corresponds to our symbol 6"). In the following, we suppress the irrelevant flavor index 

6 b ++l, |0} ^fi, c k ^c±(±k) 1 N++j±(0), 

F (R±) T \ b(p) -u ^j±(±p), P >0, (B.9) 
i>(x) <r> v / 27rV±(=Fx; e), a <H- e. 

B.2 Proofs 

The flavor index set in Section IB .11 is now identified with the set of pairs: 

X = {(s,x) : s = ±, x e Aid}- (B.10) 



Proof of Proposition \2.5[ It is easy to check that the relations in (IB. 1[) and (IB.2p are 
equivalent to the operators 

<Mk) = Cr,sMe~ ik - sX (B.ll) 

xeAiD 

and ipl s (k) = ip r s (k)t, k = k + e + + &;_e_, satisfying the relations in (12. 2p and (12. 3p . More- 
over, with this identification, the operators 

JrM= ^ !/ ' ~'' (B.12) 

xsAid 

with p = p+e + + p_e_, are identical with the nodal densities in (12. 6p . Straightforward 
computations show that Parts (a) and (b) of Proposition IB. II are equivalent to Parts (a) 
and (b) of Proposition 12. 5| respectively. □ 



20 



Note that all fermions in [3T] have chirality r — + (at least in the first part of the text). 
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Proof of Lemma \2. 7[ Using the identification in fIB.lOj) . Part (c) of Proposition IB. II is 



obviously equivalent to Part (a) of Lemma 12.71 Furthermore, the inverse of the relations in 

pmp is 

Jr,sM = 2 Yl JrAP)^ (B.13) 

p_ s eA* D 

and allows to express any operator j r>S}X (— rp s ) with p s > as a finite linear combination of 
the boson creation operators 6j(p) in (12.261) . Thus Part (d) of Proposition IB. II implies that 
the states in (I2.37P are a complete basis in J 7 . Using the canonical commutator relations 
of the boson operators b^\p), and the relations in ( 1B.4I) and ( IB. 51) . it is easy to check that 
the states in (I2.37P are orthonormal, and this implies Part (b) of Lemma 12.71 Part (c) of 
this lemma is proved in a similar manner. □ 



Proof of Proposition [Ql Use the identification in ( IB. 101) . insert ( IB. 131) in ( 1B.7[) . and then 



use llfj.sil m 



flEJU). □ 



C Solution of the Mattis model: additional details 

In this appendix, we give further details on how to obtain the results in Section |3j 



C.l Neutral bosons 

To fix our notation, and for easy reference, we collect some simple results on neutral bosons 
needed in Appendix IC.2I 

Consider the Hamiltonian 

N 

h = Y, ^Wmbm + blmb-rn + l) {^m > 0) (C.l) 

m=l 

acting on a boson Fock space defined by annihilation (creation) operators bm , labeled by 
quantum numbers m = ±1, ±2, . . . , ±iV, and a vacuum vector Qb, with 

[b m , K] = 0, [b m , b ] n ] = 8 mtW b m n B = 0. (C.2) 
Neutral boson operators P m and Z m are defined by 

/a° 1 

P m = -i\ ^(b m -bl m ), Z m = -==(& m + &Lj (C.3) 
V 2 Mo 

Y \m\ 

such that 

[^rni Pn] = i^m,n, [^mi ^n] = [-^mJ -^nl = ^' ^rn = P-mi = ^ -m- (C-4) 
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In terms of these, the Hamiltonian in ( IC.lj) becomes 



N 

h = Y, ( P l P m + (>Sn) 2 ZlZ m ) . (C.5) 



m=l 



With these preliminary definitions at hand, we introduce a more general boson Hamil- 
tonian of the form 



N N 



h = £ (PlA m , n P n + ZlB m>n Z n ) + {zlK m + KlZ m ) (C.6) 



m,n=l m=l 



with A = (A mn ) and B = (-B m) „) self-adjoint, positive definite, N x N matrices, and 
K = (K m ) a complex vector. It is useful to write this Hamiltonian in matrix form as 
follows, 

h = P f AP + Z f BZ + Z f K + K f Z. (C.7) 

We define C = A^BA 1 / 2 . Since the matrix C is self-adjoint, it can be diagonalized by 
a unitary matrix U. We write A 2 = U^CU with A the (positive) diagonal matrix whose 
elements are the eigenvalues of C 1 / 2 . 

Our aim for the remainder of this section is to explicitly construct a unitary operator 
that transforms h in flC.6j) into a diagonal form similar to the one in (1C.1I) ; see Equation 



flC.14p below. This is done in two steps stated in Lemmata 10. 1 1 and 10.21 below. 

In a first step, we construct a unitary operator W that transforms h in flC.61) into a form 
similar to the one in (IC.5I) . 

Lemma C.l. Assume that U = e M with = — M. Then 

jji __ e ,i(KtB- 1 P+ptB- 1 K) e -i(ptln(A)Z+Ztln(A)P)/2 e -i(P+MZ-ZtMP) (C 8) 

is a unitary operator implementing a canonical transformation (Z,P) — > (Z, P), i.e. P m = 
U'P m U /J< and Z m = W Z m U ,J[ , and this canonical transformation is given by 

P = A~ 1/2 UP, Z = A 1/2 UZ - B^K. (C.9) 

Moreover, this transformation applied to (C.Tj) yields 



U'^hU' = P^P + Z t A 2 Z — K^B^K. (CIO) 

(Proof: See the end of this section.) 

The diagonalization of the Hamiltonian in fIC.lOj) is now straightforward. Introducing 
boson annihilation (creation) operators bm through (cf. (10.30 ) 



P m = - b - m )' Z m = ^f^iKn + 0L m ), (C.ll) 
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with A m the eigenvalues of C 1//2 (= diagonal elements of A), one obtains 



N 



U«hU' = Kiftfim + b-J-m + 1) " K t B~ 1 K. (C.12) 



m=l 



In the second and final step, we define the unitary operator IA" that transforms between 
the boson operators bm and bm- 

Lemma C.2. Introduce the unitary operator 

/ - x A — A 

U" = exp (- J2 » m {tfJ-m ~ b- m bj), tanh(/i m ) = m " (C.13) 

m=l A m + A m 

and set U = U'U" . Then 

N 

rfW = Ki^lK + O-™ + 1) " K t B _1 K. (C.14) 



m=l 



(Proof: See the end of this section.) 

Using this lemma, one can construct all exact eigenstates and corresponding eigenvalues 
of the Hamiltonian in ( 1C.6I) . 



From the last two relations in (1C.4j) follows that ( 1C.6|) can equivalently be written as 



N N 



h — (P- m A nm P_ n + Z l _ m B nm Z_^ + [z]_ m K_ m + K^_ m Z_ m j (C.15) 



771,71=1 771=1 



with K_ m = K m etc. 



As noted in Section IC.21 the boson part of the Mattis Hamiltonian can be expressed 



as an average of the right-hand sides of ( 1C.6I) and ( 1C.15I) . One can thus diagonalize the 



boson part by a unitary operator U, using the results above. As explained at the end of 
Section I3.1[ this reduces correlation functions of the Mattis model to expectation values of 
products of simple boson operators for a system of decoupled harmonic oscillators (up to 
the zero-mode terms which require a separate treatment; see Section l"C3j) . 



Proof of Lemma \C. 11 The Hamiltonian in (IC.6I) can be written as h = P^P + Z"I~A 2 Z — 
K^B X K with 

P = U f A 1/2 P, Z = U f A- 1/2 (Z + B~ 1 K) (C.16) 

equivalent to ( 1C.9I) . 

We consider the unitary operators 

U x {a) = e l &* +pia \ W 2 (L) = e ptLZ - ztLtp (C.17) 
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for a complex iV- vector a = (a m ) and a A^ x A^-matrix L = (L m>n ). These operators 
implement canonical transformations as follows, 

Ui(a)PM(ay = P m , W 2 (L)P m W 2 (L)t = (e~ iLt P) m , 

Ux{ol)ZJAx{o^ = Z m + a m} U 2 {V)Z m U^ = (e~ iL Z) m . 

(For the bottom left relation in ( 1C.18j) consider Z m (t) = K\{ta) Z m Ui(tcx.y with a real 
parameter t. We compute dZ m (t)/dt = Ui(toc)i[a^P + P^ot, Z m ]Ui{tcx)^ = a m using (1C4[) . 
This together with Z m (0) = Z m implies Z m (t) = Z m + ta m which, for t — 1, gives the 
result. The other three relations in ( 1C.18I) can be proved in a similar manner.) 

By inspection, we see that the canonical transformation (Z,P) — > (Z,P) in (1C 16|) cor- 
responds to a sequence of three canonical transformations implemented by the unitary op- 
erators W 2 (-iM), W 2 (-iln(A)/2), and ^(B^K); see flUT7j) - (IUT8|) . Thus, the full trans- 
formation (1C.16I) is implemented by the operator W = Wi(B _1 K)W 2 (-iln(A)/2)W 2 (-iM), 
equal to ffClSjl . □ 

Proof of Lemma I CM Computations similar to those in the proof of Lemma IC.ll give 

U"b m U" ] = cosh(/X| m |)6 m + sinh(//| m |)6L m . (C.19) 
Combining ( \C3\i with fIC.lip . and using the fact that 





cosh(^ |m| ) = 1T T + ^ sinh(^ H ) = , - . , (C.20) 



one sees that the right-hand side of (IC.19P is equal to b m . Equation (IC.14p follows from 



C.l.l Correlation functions 



We give some well-known results on (free) boson correlation functions needed in Section [3j 

Let be boson operators satisfying (1C.2I) and Hb = Y2m u m^ln^mi unlike before, in 
this subsection Yl m means t ne sum over m = il> =t2, . . . , ±iV, and similarly for products 
Y\ m - Denote by (A)b,/3 the usual thermal expectation value of a boson operator A with 
respect to that free boson Hamiltonian, i.e. (A)b,p = Ttb(g~ i3Hb A) jZ B ,p with Tr B the 
trace over the boson Fock space and 

m 

see e.g. Appendix C in [21] for more details. 
Lemma C.3. (a) For j = 1,2, let 

^3 = / j {J-hmPm + Jj,mPm) (C.22) 

III 
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with complex numbers J® and Jj m - Then 

{J\J2)b,P = J[A + Yl { J l,m J 2,m x _ e _^ m + Jtm-Km _ J ■ ( C ' 23 ) 
m 

(b) For j = 1,2,...,N, let 

K t = J2 K trJ>L Kr=J2 tq m b m (C.24) 

m m 

with complex numbers K im . Then 




(C.25) 



Outline of proof. Part (a) is a simple consequence of the following well-known results, 
(K b L) B,p = S n , m /{1 ~ e~ PWm ) and (b n b m ) B ,p = (b m ) B ,p = 0. Part (b) can be proved using 
well-known identities collected and proved in [2T], Appendix C, for example. To be specific, 
writing Kj = + Kj , one finds 

e *i e K2 ... e KN = e E i K * e^ ^' K ^ = 

using Theorem 2, Equation (C4) in [21], and (1C.2j) . respectively. The result follows from 
the relation 

( e £^) = e ^ jth (W/2 (a27) 
(see Theorem 4, Equation (CIO) in |21j). applying (|C.23p . and by collecting terms. □ 



Remark C.4. In the proof above we suppressed technical details concerning the unbound- 
edness of the operators Kp. It is well-known how to take care of these details using Weyl 
operators; see e.g. [31], Proposition 2.1 ff. 



C.2 Proofs 

We give the proofs of Theorem 13.11 and Proposition 13.41 stated in Section [3j 

The Mattis Hamiltonian H3 .31) can be written as a sum of Fourier modes H(p) that 
mutually commute. In order to diagonalize the Hamiltonian, one can either treat each 
mode H(p) separately using the results in Appendix IC.ll for N = 2 or N = 1 (the latter 
for the cases p = p±e±, the former otherwise), or treat all the modes in one go. Although 
we have opted to present the latter approach here for conciseness, the careful reader might 
prefer to write out the explicit details of the former approach. Note that the matrix B 
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in Appendix IC.ll is non-trivial only in the cases when '"f2X{p)p+P- 0- Moreover, K is 
non-zero only if p = p±e±. 



Proof of Theorem\3Jl From fl2T26|) and follows that 



fi s (p) = -i^y ^(6.(P) -6J(-P)), *.(p) = 2^WP) + ^(-P)) ( C - 28 ) 
analogous to (1C.3jl . 

The Hamiltonian in fl3.3p -( l3~5l) is of the form H = vp :h: +Hq^ with h the average of 
the Hamiltonians in (1C.6j) and (1C.15j) if we identify m with (s, p), s = ±, p G A*, and set 

27T - 27T « 27T - 

P fl , p = — n a (p), Z S:P = j-$ s (p), K s , p = —E s (p) (C.29) 

with — m corresponding to (s, — p). The matrices A and B have block diagonal form, i.e. 
A s> p, s ',p' — ^4s,s'(p)^p,p' an d similarly for B, with 

As,s>(p) = (1 - 7iX(p))^,s', B s>s/ (p) = ((1 + x(p)7i)<W + x(p)72*«,-*')p*Pa'- (C.30) 

It follows that the matrix elements of U in Lemma IC. II also have the form U SiS i(p)8 Pj pi, and 
similarly for M. We note that A and B commute, and thus C = AB (see Section \C.l\i . 
When 72X(p)p+P- 7^ 0) the diagonalization of C reduces to that of 2 x 2 matrices. When 
12X(p)p+P- — 0, the block B s s /(p) of the matrix B has only diagonal components s = 
s' = ±, as is the case for all other matrices A, C, U, and M. Furthermore, we can treat 
K as a complex vector since its components commute with those of P and Z. With the 
identification in (]C.29|) . the commutator relations in (13 .2p are equivalent to the ones in 

By straightforward computations we find that the eigenvalues X s , p of C 1 / 2 = (AB) 1 / 2 , 

multiplied by Vp, are equal to w s (p) in (13. 7p . and the result for U s>s i(p) is given in (|3.25|) . 
Moreover, 

(B^K)^ = ^ 5s( P } - = ~-^± X (p)Q + MPs)5 P - s ,o (C31) 
L (1 +7iX(p))PS L 1 + liPs 

for p G A*, and 

K ^ lK = YTvlH £ {^) 2 x(p*-s)Q + A-p)Q+Ap), (c.32) 

s=± P e^z\{o} 

and for the unitary operator in ( 1C.8I) we obtain (see also Footnote UM on page fT6i) 

w'=ex P (i;)r;r;(^) 2 (^ 

\ s=± pe A; 



X 



ex P( -*£ E (^) a in(l-7iX(p))(fll(p)^(p) +^(p)n s (p))/4 ) (C.33) 



: P eA* 



ex P h'E E (y) 2M .- s (p)(nl(p)^(p) - $l(p)n- s (p))/2 

V s=± pe A* / 
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with 

M s ,Ap) = -s6 a> .s> arctan([/_+(p)/?7 ++ (p)) (C.34) 
and U s>s ,{p) in (jX25]l . Thus (lOTOD and flUT2|) give 



W*:/i:W = E E (^K(P)6.(P) + 1/2] - \p.\/2) - K t B~ 1 K (C.35) 
where we write b^\p) = &i!p; we used u s (p) = u s (— p) and 



4<p) = -^^(to>-«H»)) 

^ p »4iSf (i ' (p)+ ' :| - pl) 



(C.36) 



implied by (IC.llj) . The additional constants —\p s \/2 in ( 1C.35P are due to the normal 
ordering. 

Following Lemma IC.2t we set 

-27T\2 



U" = exp ( i E (t) ^ (P) (^^^(P) + *!(P)n.(p))/2 ) (C.37) 

with 



P6A* 



tanhfe(p)) = "'' P |:" r ' P '| (C.38) 

and where we have used (10.281) . Finally, we define the unitary operator that diagonalizes 
the Mattis Hamiltonian as the product of ( 10.33[) and (10.37|) . 

U = U'U", (C.39) 

such that 

uU - h '- U = E E (^K(pMp) + 1/2] - N/2) - K t B~ 1 K. (C.40) 

Since H = v F : h : +^ 0) , (HT40D implies (ESDif with ifg = - w F K t B _1 K and £ 
in (IBTTTTI) . U sing (g 5) and ffOT32]) . together with the Parseval relation (12T39D . (PX42D . and 

4t,,(0) = ^a/(27r) E^eAu, <2±, s 0), one obtains (O- □ 



Proof of Proposition 3.4 Since N r>s (x) obviously commutes with W in (1C.39I) and H in 
(BID , (BUBD is trivial. 



For the Klein factor, note that U'^R r , s {x)U' = e- iD R TiS (x) e iD (recall W in (10331 ) with 

D = E E (f) 2 2 fi+\ b 2 ^ I(p)fi - (p) + ft * (p) ^ (p) ) (c - 41) 



s : peA* 
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and H a (p) in (13. jj) . To compute this we define f(z) = e lzD R rtS (x) e lzD , z e C, and find 
df(z)/dz using 

[Q+ ia (p),i2r,«(a;)] = ry/K/nRr^x) e~ lpx /2 (the latter follows from fTOUj) . 
(I2.32p . (I2.34p . and (I2.42p ). The resulting differential equation for f(z) can be solved with 
the condition /(0) = R rtS (x), and setting z = 1 gives 



U*Rr >s {x)U = ^(x) exp (-r 2 ^ ^| £ (t) 2 ^-^-^ 6 "^)- (C - 42) 
Using (1C.28j) and (1C.36I) . one sees that 



U"^ s {p)U" = J^rfl s (p), U"%(p)U" = J^4&s(p), (C.43) 
V v F \ps\ y w s (p) 

from which R T ^ s (x)U then follows. In a similar manner we compute 



e [Ht 



R r ,s( x ) ^ = e ikrAx)t/2 R r ,s( x ) e ik ^ x)t/2 (C.44) 



with k r>s (x) in fl5720|> . We recall fT0728]l and obtain (|3TT5]l - (^2Tp by straightforward com- 
putations. 

Lemma EH f!U^9|) . and flU3B imply 

Wn.(p)W =(i - 7 ix(p)r 1/2 £ ^,,'(p)n S '( P ) 

s'=± 

U^ a {p)W =(1 - 7iX(p)) V2 X) ^(P)^-'(P) (C.45) 



s'=± 



+ ; ^ X(p)— Q+,-s(p s )^_ s ,o 
1 + 7i Ps 



for s = ±. Using (1C.43I) . J r ,s(p) = A/7r/a(ip s ( l ) s (p) — rfl s (p)) for p s ^ (this follows from 

(EED), (PSJ, (J325K, and 

e i# *6 s (p) e- i# * = 6 s (p) e - iaJs(p)t , (C.46) 
we obtain the result in (l3T23|) - fl3~24l ). □ 



C.3 Zero modes 

We show that the contribution from zero modes to the free energy density and any corre- 
lation function is irrelevant in the IR limit. Much of what is stated in this section has its 
origin in the following result. 

Lemma C.5. Let H = H T be a real positive definite M x M matrix, m G R M , an$3 

Z= e' uTHu+imTu . (C.47) 

21 Note that Z is (a variant of) the Riemann Theta function, albeit using an unconventional notation. 
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Assume 

\ d J min xTH " lx > A ^ ^ max I (H -1 m) I < 2ttA. (C.48) 



T/ien 

Z 



J 

with 



1 < - < 1 + (MA 2 e- A " 2 ) (C.49) 



Note that H _1 exists since H is positive definite, and that 1/A is equal to the largest 
eigenvalue of H. Furthermore, the exponential damping in (IC.49I) shows that (IC.47P can 
be estimated, to leading order in A, by replacing the sum by a Gaussian integral. 

Outline of proof . Using the Poisson summation formula and ( 1C.50I) . the sum (1C.47I) can be 
rewritten as 



71 '_ ST^ e -|(m+27rn) T H- 1 (m+27rn) _ j g — 7r 2 n' r H- 1 n-7rm T H- 1 n /q ^\ 

det(H) ^ ^ ' 

The last sum is bounded by 

1 < e- 2nTH " ln -™ TH " n < (j2 e'^+^A (C.52) 

with A and /z in (IC48[) . For /i < 2ttX, it holds that 



oo 

2, 



j2 e-A^n'+Hnl < 1 + 2 e - Xw2+ ^ + e^ 2/(4A) e _W (C.53) 



nSZ \ n=l 

and, if A > 2/n 2 , one finds (see for example [19], p. 153-156) 



oo 

^ e - w = C>(A 2 e- A7r2 ) (C.54) 



n=l 



from which flC.491) straightforwardly follows. □ 



C.3.1 Free energy density 

We derive the right-hand side of (I3.14p . The zero-mode contribution to the partition func- 
tion is (cf. ( El ) 

Zq,P = J>°>- ^No,,} (C55) 
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with Hq in (13. 9p . This can be written in the form (IC.47j) with m = and M = AL/a. 
From (13. 9p . one sees that A = O(L), and it is immediately clear from Lemma [C. 5 1 that the 
leading order contribution to ( 1C.55|) is obtained by the corresponding Gaussian integral. A 
straightforward calculation yields (using f]C49|) ) 

(j \ 2L/a 
^—Jj) (I + O (L 3 e~ cL )) (C.56) 

with A and vf in (13.81) . and c > some constant independent of L. This also implies that 
the IR limit of the zero-mode contribution to the free energy density vanishes 

^f = -JT2 H^QJi) = ln L ) ■ ( C - 57 ) 
C.3.2 Correlation functions 

The irrelevance of contributions from zero modes to thermal expectation values in the IR 
limit is a simple consequence of the following: 

Lemma C.6. Let m r s (x) G M, with r, s = ± and x G A^. Then 

T £ £ mr,a(x)N T ,s(x) \ 

/ Q,P 

-s^r E £ E ( 3I 4 J + I ^) mP , s W^, s W (C.58) 

x e ^' =± s=± "■"' 6A "> 1 (1 + O (L 3 e- cL ) ) 

wzt/i c > some constant. 

Proof. By definition, 

i; E E m r , s (a;)Af r , s (x)\ 1 ^-^ -P^Q + T E E »nr,»(*K, t (i!) 

r, S =±x£A 1£) \ = \ g V=±rtA 1D _ (C.59) 



As in Section IG.3.1[ we can apply Lemma IG.5I The result follows by replacing sums by 
Gaussian integrals. □ 

All iV-point correlation functions involving zero modes can be generated from ( 1C.58I) 
by differentiation with respect to the m rjS (x). We have, for instance, the important special 

case 



1 1 / / 1 T\T2 \ 



a 72 



f<U,«) (l + 0(L 3 e" ci )) 



(C.60) 



LA(l + ll ) 2 

whose IR limit clearly vanishes. We also note from the above lemma 

/ r £ E m r,s(x)N r , 3 (x)\ 
lim ( e '•■•=±-eAi D \ = i (c.61) 

L^oo \ / Q,f3 

if all but a few of the m r s (x) are non-zero (i.e. the number of non-zero terms is 0(L°)). 
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C.3.3 Computation details 



We end this section by verifying the claims made regarding zero-mode contributions in 
connection with Results 13.51 and 13. 6[ respectively. 



On Result HO . From the scaling with L in (13. 2 7ft and comparing with ( 1C.60I) . it is clear 



that ( J° 1Sl (xi, t%; e) </° 2)S2 (x2, e)y = 0{L x ). Thus, the zero modes do not contribute 

to the density- density correlation functions in the IR limit. The same holds true for higher- 
order density correlation functions, as seen by successive differentiation of the generating 
function f l058|) . □ 

On Result \3.6l We want to compute the zero-mode expectation value (13. 36ft for the oper- 
ators S q s (x,t; e) in (I3.32p . One can move all terms of the form e~ iqKr ' s( - :x '' t ' t ^ 2 to the right 
by successively applying 

e- iqK °»^' t;e)/2 R r , iS ,{y) r ' q ' = R r ,, s ,{y) r ' q ' e ^<i K U^)/ 2 e m'^ r , a y,A^v^)/L ( C 62 ) 

where 

£ r ,s,r>, s >(x, y, t; e) = ^- ( [(1 + ^ A + rr'{l - 71)] 8 x _ a>y 8 S]Sl + y7 2 

r 72 A (C - 63) 

x [-7—. S a ,a' + S s -s'\ ) - rx s S r yS s ^S x _ a>y + r5 s _ s ,f(x s - y; e) 

1 + 71 / 

(this follows from (13.331) and the commutation relations in ( 12. 34ft ). Then (I3.36f) becomes 

X q 1 T Ej=i ^ r j> s j> r j> s j ( x j' x j,-sj i*i; <E )+ 2i f' J2i<j<k<N 1i1^ r r a r r k- a k ( x j' x k > -s k ,tj;e) 



The second factor becomes 1 in the IR limit (cf. ( 1C.61I) and ( 13. 33ft ). The same holds true 
for the third factor, leaving us with (13. 37ft . □ 



D Quantum field theory limit: additional details 

D.l Free energy 

We outline the computations leading to Result 14.11 The zero modes do not contribute in 
the IR limit by ( IC.57I) . and it suffices to only consider the boson part given in ( I3.13p . We 
separate the latter into terms with quantitatively different momentum dependence for the 
dispersion in (13. 7p . writing VLb,p = + with (in the IR limit) 

TV /a TV fa 

L->n h = p.v. J J If^jMi - .-*«). 

-7r/0— 7r/o ' /j-j j\ 

OO \ " / 

lim L~ 2 n> g =4- / ^Iln(l- e"^), 

7r/5 
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where the Cauchy principal value for the improper double integral is defined as 



tx la 



P.V. 



dp A 



lim 

L— >oo 



-7r/a 



-2-k/L 


7r /a 


J d P± (-- 


'>+/ 


-7r/a 


2-k/L 



(D.2) 



The second term in (ID.lj) is exponentially suppressed 0(e~ l3vF / a ) and thus vanishes in the 
UV limit (although see the footnote on page [24]) . For the first term, we consider initially 
the simpler case of 72 = 0. By inserting (13. 7j) . using the symmetry of the integrand, and 
rescaling the integration variable, we get 



2n< 



lim L~ Z Q 

L— >oo 



B,f3 



TTVFa/3 2 



(72 = 0) 



(D.3) 



which, in the (properly scaled) UV limit, becomes 



2n< 



lim lim aL £1 



7T 



B,/3 



36 F( 5 2 



(72 = 0), 



using that 



00 

/ x z In (1 - e~ x ) dx = -T(z + l)((z + 2) (&(z) > -1) 



(D.4) 



(D.5) 



and £(2) = tt 2 /6. For the computation with non-zero 72, it is useful to introduce the 
effective dispersion in (14. 2 p and a corresponding free energy contribution tig*. The QFT 

limit of tig a is computed as in ( 1D.3j) ( IT3T41) . and this gives the right-hand side of (14. ip . 
In order to prove Result 14.11 it is enough to show that the remainder term Qp a — Q% a is 
sub leading in a and therefore does not contribute to (14. ip . In the IR limit, 



k/cl k/cl 



L— >oo 



lim L~ 2 



B,f3 



L B,/3 



P.V. 



-TV /a —rc/a 



d 2 P 1 

(2^ 



_ p -f3u 3 (p) 



(D.6) 



Symmetry of the integrand allows us to restrict the integration domain to the following 
triangular region: < p + < tx /a and p_ < p + . Changing to polar coordinates gives 



L~ >oo 



lim L- 2 Qi a -tti 



7r/4 7r/(acos(6>)) 



with 



lim 

(2?r) e^0+ 



w±(p) = v F pg±(9), g±(6) 







-(3w s (p) 



-/3S s (p) 



(D.7) 



1 ± Wl - A sin (20) 



(D.8) 



and similarly for u)±(p). We note that both the numerator and the denominator in the 
argument of the logarithm in f 1D.7|) become zero as p — > (for s = ±) or 9 — > (for s = — ). 
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However, their ratio is non-zero and finite in these limits (this is the reason for the specific 
choice of (14.2j) ). and the double integral therefore has a well-defined limit as e — > + . 



To finish the proof of Result 14.14 we show that the right-hand side of (ID . 7[) remains finite 
as a — > + . The expression ( ID.7|) is of the form J f(8)d9, with / given as an integral 
over p. By the mean value theorem, there is a 6* in (0, vr/4) (possibly depending on a) such 
that /; 7 V(0)d0 = / (0*)tt/4, i.e. 



lim L Q 



L— >oo 



7r/(a cos(#*)) , 

^ s_± 



-/8w«(p) 



pdp, (D.9) 



which, after changing integration variables, can be written as 

0VF9s{O*)n/ (a cos(0*)) 



Hm l- 2 fn^-n^ 



1 1 1 



2tt S 2 ^3 



1 



dx — 



1 



8=± ' 

PvfVA tan(0* )7r/a 



dx 



Acos 2 (#*) ./ Asm 2 

If 6* /a remains finite in the UV limit, it follows that 
lim lim L~ 2 ( tti 



dx 



x 



In (1 - e- x ) 



a->-0+ L->-oo 



1 B,/3 



where we used (1D.5j) to compute the integrals. Otherwise, one finds 

1 1 



lim lim L~ 2 VL< 

a->-0+ L^oo v ■'- 



1 

2ir 



X 



lim [ 

a->o+ l£± 1 + s^l - A sin 2 (26*) Acos 



Asm 2 



(D.10) 



(D.ll) 



(D.12) 



Using trigonometric identities, one shows that the terms in parenthesis on the right-hand 
side above add identically to zero for all 9*. This finishes the proof that the remainder is 
subleading in a. 



D.2 Fermion two-point functions 



Result 13.61 and Equation (I3.4ip imply that the IR limit of the fermion two-point functions 
at zero-temperature are 



lim (< s (x,t)e, s ,(0,0) 



x lim hm — — r— 

l-s>oo £-s>o+ 27rae G rjSiriS (0, 0; e) 



(D.13) 
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with (cf. flSSH ) 



peA* 



s'=± 



(D.14) 



D.2.1 Proof of Result H31 



For 72 = and |p±| < vr/a the formulas in (13. Tj) and (I3.24p simplify to 

w a (p) = ufIp.I, <' s (p) = iy^<5 s , s /^jp7| (sgn(p s )/B + r/V#) (D.15) 
with vp in ( 13. 8 p and B = v F (l — ji)/v F as in ( 14. 7p . Inserting this into (ID. 14^ gives 



In G r . o [XcQc , t, 6j 



|p_ s |<7r/a\ p s ^0 '^ S| 



x e -WF\Ps\t Q ipsx a e -e\p a \ 



|p s |>7r/a 



(if + rsgn(p s )) e 



with if in ( 14 .TP and p± G (27r/L)Z; the second p s -sum in ( 1D.16h is to account for the fact 
that ( 1D.15P holds true only for |p s | < ir/a, and u s (p) = Vp\p s \ and |v^ s (p)| 2 = 5 SiS /|p s |(l + 
rsgn(p s ))/(47r) for \p s \ > ir/a. The sum over p_ s in (ID. 16[) cancels the overall factor a/L. 
The first sum over p s can be computed using the series Y^n=i e~ zn / n — — ln(l — e~ z ) = 
-\n(z) + O(z) for $l(z) > 0, leading to 



In 



2tt 



K 



(e - \{rx s - v F t))({e + w F tf + (x s f) 



2\(^-l)/2 



(D.17) 



up to 0(L x ) corrections, while the second sum over p s can be approximated, up to cor- 
rections of 0(L _1 ), by the integral 



•'f ^ r'=± 



-p^+rCpi— ir'ra; s ) 



(?(e + rM - irx s )) V" (if + r')#i (?(e + iv F t - ir'rx s )) 



r'=± 



(D.18) 



with defined in (14. 4p . Putting all this together, we arrive at 



A 



1 



V27ry e - i(rx s - v F t) \(e + iv F t) 2 + (x s ) 



1 



(K-l)/2 



a 



(l(e-i(rx s -v F t))) 



(K+l)/2 



a 



(l(e + i(rx s +v F t))) 



(A'-l)/2 



(D.19) 



a(l(e-i(rx s -v F t))) 
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{1 + OiL- 1 )) 



with a(z) defined in (14.41) . The function E\ has the asymptotic behavior (see e.g. [50] , 
Equations 5.1.11 and 5.1.52) 

E 1 (z) = - 7 - ln(*) + 0(z), E 1 (z) = (l + O ((z + l)- 2 )) , (D.20) 
implying (14. 5p . It follows that, for 6<(i, 

G r , w (0, 0; e) = (J£j (e 7 £)*~^ (l + O (e/a) + 0(L- 1 )) . (D.21) 



Inserting this and flD.19j) into flD.13j) . we obtain (I4.6p -( H~8 



D.2.2 Proof of Result l4~4l 

The computation of (ID. 13[) (ID. 14[) for non-zero 72 mimics the computation of the free 
energy in Appendix ID. II We add and subtract to the non-trivial part (i.e. In Gf s r J of 
( 1D~T4|) a term 

l*C%„(x.e. t t-,e)*a £ (|) " ^ 1 + rp ' e^ e^^ e^' (D.22) 



PGA* 
|p s |<7r/a 



with w s (p) in (14 .2p and in (14. lip . We note that this expression follows from that of 
In Gf s r s by everywhere replacing the original dispersion relation uj± (p) with the effective 
dispersion u)±(p) (i.e. making the replacement in the coefficients f!3.24j) — (13. 25j) for the case 
72 = 0, and in the time-dependent part of (1D.14P ). Result 14.41 can then be obtained by 



analyzing the QFT limits of G r ^ rjS = G< srs G> s ^ s and E r ^ rjS = Gf JGf separately. 
The former can be computed as in Appendix ID. 2 11 and one finds expressions identical to 
(1D.19P and (ID.2ip . except with v F replaced by Vp\/A and K now given in (14. lip . In the 
regime (14. 3p . we obtain 



2^/ e - \{rx s - v F \fAt) 

(K-l)/2 

(1 + 0(S))(1 + 0(L- 1 )) 



(D.23) 



{e + w F VAt) 2 + (x s )\ 
using (the equivalent of) (ID. 19[) and the second asymptotic relation in (I4.5p . 
The error term in the IR limit is given by 



TV /a 7r I a 



Ymv^ In E rtS ^ s (x s e s , t; e) =P.V. J J d 2 p a ( ^ 



|2 



\ V r,8\P)\ c -iov(p)f 



\Pi 



1 2 



-n/a-./a Vs ' =± (D.24) 

l — (K \p,\ + rp s ) e-^-W* j e ip ° Xs e" e|Ps| 

\Ps\ J 



47T \p. 
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with the Cauchy principal value defined in ( 1D.2j) . We will see that, in the limit a — > + , 
this error term yields a finite constant (depending only on the coupling parameters) if 
x s = t = 0, while it vanishes in the regime (14. 3p . Transforming to polar coordinates and 
using the independence of E rja ,r,s(0, 0; e) on r and s, we find by integrating over the radial 
coordinate and using the rotational symmetry of the integrand (we set, without loss of 
generality, s = + on the right-hand side of (1D.24|) ) 



lim \n E r s . r s ( 0,0:0 ) = lim 

L-s-oo ' ' ' v e ->0+ 



tt/4 tt/2-e 

de r do 



cos(#) J sin(6 l ) 

tt/4 



COS 2 (9) 



(D.25) 



1 ^ 4 I Vl--4sin ! (29) / V »«'(«) B ) 1 



vS'=± 



with g s '(6) in (ID.8p . The limit e — >■ + of the right-hand side exists (the integrand remains 
bounded as 6 — )■ + and 6 — > (n/2)~) and it is conveniently denoted by — In (C(7i, 72)), he. 

C(7i,7 2 ) = expf- lim In £ w , s (0, 0; 0+)) . (D.26) 

We observe that C(7i,0) = 1. 

We proceed by showing that (1D.24j) vanishes in the regime (14. 3 p in the limit a — > + . 
We only discuss the case s = — (since the case s = + is nearly identical). Furthermore, it is 
enough to restrict the integration domain on the right-hand side of (ID.24j) to the following 
triangular region, < p + < 71 /d and p_ < p + (the other seven triangles are analyzed in the 
same way). Transforming to polar coordinates and integrating over the radial coordinate, 
we find 

vr/4 

f ,„ a I 1 / s'cos(20) \ fBsm 2 (6) g s 

lim / d9 1 ' 



e i7r(sin(6»)a; s -'D F g s ;(e)t)/(acos(6')) _ ^ e i7r tan(0) (a; s — •up </~At)/a _ j\ 
^ _|_ ^ J 

sin(0)z s - t5.F3 s / (0)i ' x s -vpy/At / 

(D.27) 

The integrand remains bounded in the whole integration domain as e — > + . We can 
therefore apply the mean value theorem (cf. Appendix (ID.ip ) to its real and imaginary 
parts separately. The end result is d multiplied by a bounded complex number (the bound 
is independent of a). It follows that the UV limit of the above expression vanishes, and, 
more generally 

lim lim \nE rjS>r)S (x s e s , t; e) — 0. (D.28) 

a— >0+ L— >oo 



Inserting all results above into 
e^nL \K-i 

a->0+ L^oo 



lim lim — - — s (x,t)i> g r , s ,(0,0)) = 5 q - q >5 r y5 StS > 

L^oo V a / \ / 00 

1. 1. 1. r- / G^TtLq K— 1 1 G r s r s (x s G s ,t,6)E r s r s (x s G s ,t,€) 

x lim lim lim o x _ s ,o( — z — ) 



(D.29) 



a^ + l->oo e -^o+ s ' a 2/Tae G> s rv ,(0, 0; e) £ r , s , ns (0, 0; e 

and using that S x _ 3 o/d ~^ 8{x-s) as a — )• + yields Result 14.41 
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D.3 Density-density correlation function for 72 = 



Inserting (1D.15|) in Result I3.5[ we find after a trivial integration and by adding and sub- 



tracting terms 

lim a(j rs (x,t;e)J r > s /(0,0;e)) = 5 ss >-5 x _ sfi — — 

L->oo 00 a 4(27T 



OO 



x ( J dp \p\ (B + rr'B- 1 + sgn(p) (r + r')) q-^fW e -e\v\ jP*s + J dp \ p \ ^ ^ 

-00 \p\>n/a 

x (<5 r , r ,2 (1 + sgn(p)r) e^^* - (B + rr'B- 1 + sgn(p) (r + r')) e - ie 'lfl*) q^-^ 

Computing the first integral above gives the first two terms within the parenthesis in ( I4.14p . 
For the second integral, we scale integration variables and use 



00 

ai ( z ) = Jte~ zt dt = + 0) (D.31) 

1 



(see [50], Equations 5.1.5 and 5.1.8). This gives the last term within the parenthesis in 

dm. 



E Correction terms 



E.l Derivation 

We consider the tight binding band relations e(k) in (1A.3|) and the nodal points Q rs in 
(IA.5p . We Taylor expand 

3 3 

e(Q r , s + k) = e(Q r , s ) + v F ^ ^ a j - kik a i - 1 (rk.) i - k (k-.) k + 0(S 4 |k| 4 ) (E.l) 

j=l k=0 

with vp, given in Equation (6) in [3], defined such that a^o = 1. By straightforward 
computations we find the following other non-zero coefficients, 

(1 - k) [t cos(Q) + 2t' cos(2Q)] (1 - k) [t cos(Q) + 2t'] 

^2,0 — , ■ 77\\u 1 o-t/ TTvTi ' a °' 2 



4sin(Q)[t + 2f cos(Q)] ' u ' z 4sin(Q)[t + 2t'cos(Q)] 

(1 -/t) 2 [t + 8t / cos(g)] (1 - /t) 2 t 1 ' j 

a3,0 ~ 24[t + 2t'cos(Q)] ' ai,2 ~ 8[t + 2fcos(Q)] 

with parameters t,t',Q and k defined in [3] (we used (11) in |3J). The first non-trivial term 
VFrk s on the right-hand side in ( IE. II) corresponds to ifo in (EH]). The terms ^(a^o^s + 
a3fid, 2 rkg) correspond to H 2 in ( 15. ip . Finally, the terms (00,2^ + Oii^a 2 rk s )k 2 _ s correspond 
to if 3 in (15. 2P ; we inserted 

^ A; 2 s = 2-2 cos(A;_ s a) + 0(a 4 kt s ) (E.3) 
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and ignored the 0(a 4 ArJ-terms to obtain an expression local in the "position space" defined 
in (12. lip (this corresponds to using a lattice second derivative in the x_ s -direction). 

Note that ao,2 = at the special parameter value Q = arccos(— 2t'/t), as mentioned in 
Section [5j Remark 5. 



E.2 Bosonization 



We introduce the notation * ■ ■ • x for normal ordering of general boson multilinears, defined 
as usual: 



«Jn,Ai(Pl) • ■ -JtmAn^Pn) x = j Ml) ,A CT(1) (P<x(l)) • • ■ Jr a{N) ,A a(N) (Pa(N)) (E.4) 

with the permutation a of N objects such that r a n\p a (rj < r a M)Pa(2) < • ■ • < r a{N)Va(N)- 
Note that, for boson bilinears, this normal ordering is equivalent to the one in (12.41) . 



The following is a straightforward generalization of Proposition IB . 1 ( b) . 
Proposition E.l. For & e M, 



2iL sin I — — I 

\ L J 



-i 



x 



exp(r27ri^— — (d x ) m JrA 



m=0 



(E.5) 



with V ry A.(x; e) in (B. 7) and 



(Proof: See e.g. [26j.) 



(E.6) 



A formal expansion in the parameter b of the left- and right-hand sides of (1E.5|) . together 
with an integration of both sides with respect to x, give 



L/2 L/2 

limV- / dx:V rA (x;ey(d x ) n V r , A (x;e): = / dx 
^o+ f-^ n! J J 



X . X 

x Jr,A x 



ra=0 



L/2 



-L/2 



-L/2 



L/2 
b 2 f. 1 



+ b j dx- xr2iri(j rA ) x + — / dx- x I (r27ri) (j r>A ) + ) j r ,A j 



7T 



-L/2 



-L/2 



L/2 

+ S / dx) x ( (r27ri) 3 (j r , A ) 4 - r2m(d x j rA ) 2 + 2( | )~r27ri(j r , A ) 2 ) 3 + O(ft') 



3! J ~~4 

-L/2 



(E.7) 
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where we suppress the argument in j rr A.(x). Using that, we can bosonize the correction 
term in (15. ip and obtain 



r,s=± 



H 2 = y~] vf I d 2 x x ( aa 2 , -(vra) 2 (J, 



3\L 



(E.8) 
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where we use the notation in (12. 17ft . We insert ( \1.6h and obtain ( 15. 3ft . 

Remark E.2. The bosonization of higher derivative terms has been discussed in the literature 
before; see e.g. [51] and references therein. 



F Index sets 



A s = 

a: = 

k = 
a: = 

Aid 



m L L 

x : x s G K, x_ s G aZ, < x+ < — 

' 2 ~ 2 







t( z+ 


-Y 




2/ 




7T 






T z : 


a 



7T 

a 



7T 



We collect for easy reference the definitions of various index sets used throughout the text: 

dns]) 

fl2~27D 

fl23B} 

d23I]) 
dS2Z]) 



a: 



ID 



p G A: : p s ^ 



x G aL : — — < x < — 

271 IT 7T 

p G — Z : — — < p < — 

L a a 



A 1D = <peA* m :p^0 
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